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We describe an optical apparatus designed and built to extend con- 
ventional light-scattering measurements to the very-small-angle regime. 
The present instrument covers the angular range 0.008° ^ 6 ^ 0.15° 
with an instrumental resolution (hwhm) of 0.00045° (1.6 arc seconds), 
and exhibits an exceptionally low stray-light background. The theoretical 
and practical considerations important in achieving this performance are 
analyzed in detail. Besides its primary purpose of studying long-wave- 
length (0.01 cm to 1 cm) thermally driven fluctuations, the present type of 
apparatus should also prove quite useful in other areas where long-wave- 
length perturbations must be probed, such as, (i) holographic and optical 
memory imaging, (ii) surface roughness testing, and (Hi) index of 
refraction profiling. 

I. INTRODUCTION 

Laser light scattering has, over the past decade, been developed 1-3 
into an extremely powerful tool for probing the long-wavelength 
(A/ « 2 X 10 -5 cm to 2 X 10~ 3 cm) elementary excitations of liquids, 
gases, and solids. Combined with diffraction grating, Fabry-Perot, or 
optical mixing spectrometers the technique is capable of spanning an 
impressive range of more than 13 decades in energy or frequency 
measurement. Yet, in contrast to what has happened in the field of 
inelastic X-ray scattering, 4 - 5 very little has been done to utilize 
very-small-angle (vsa) light scattering to probe longer-wavelength 
(A/ -£3 10~ 8 cm to 1 cm) excitations. With a few notable exceptions, 6-11 
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most light-scattering experiments have been limited to the scattering- 
angle range > 1°. 

There have been a number of reasons for this apparent lack of 
progress in the very-small-angle scattering regime. On the one hand, 
experimentalists in the field, encountering a seemingly divergent stray- 
light level at small angles, have assumed that attempts to work in 
the vsa region would present insurmountable problems. On the other 
hand, there did not appear to be any physical phenomena where the 
important elementary excitations were confined to the corresponding 
longer-wavelength regime. Or, in cases where they were, it seemed that 
the use of more conventional macroscopic experimental techniques 
represented a satisfactory experimental approach. 

Recently, however, there has been a resurgent interest in problems 
involving general hydrodynamic instabilities 12 both in normal liquids 
and liquid crystals. 13-26 The "critical wavelengths" involved in the 
onset of these instabilities are, in general, controlled by some macro- 
scopic dimension of the sample chamber and tend to fall in the range 
100 nm < A c < 1 cm. Light scattering is the only technique offering 
the possibility of probing these wavelengths without physically disturb- 
ing the sample and with a sensitivity sufficient to detect the thermally 
driven critical fluctuations. However, probing the excitation wave- 
length region 100 Mm < A < 1 cm requires the capability of resolving and 
detecting the scattered light at very small angles, 0.3° ^ 6 ^ 0.003°. 

This paper describes the experimental progress which has been made 
in extending the light-scattering technique to this very-small-angle, 
long-wavelength regime. 

In Section II, we describe the physical configuration of a light- 
scattering apparatus that has been constructed for use in the vsa 
region. This section also summarizes the measured performance 
characteristics of the instrument in terms of angular resolution and 
stray light. Section III is a detailed presentation of the basic diffraction 
and aberration considerations that influence the design of a vsa 
light-scattering apparatus. Section IV outlines various empirical obser- 
vations made during the course of construction of the present instru- 
ment, relating to the stray-light behavior of optical components at 
small angles. 

II. AN APPARATUS FOR VERY-SMALL-ANGLE LIGHT SCATTERING 

2.1 Introduction 

In this section, we present a general description of the physical 
layout and performance of a light-scattering apparatus that has been 
constructed for the vsa regime. The theoretical background and 
practical considerations necessary to analyze the detailed charac- 
teristics of the instrument are deferred to Sections III and IV. Although 
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designed specifically for the study of the B6nard convective instability, 
this apparatus embodies solutions to most of the problems to be 
encountered in the general small-angle light-scattering experiment. 

2.2 Performance goals 

The following performance goals were established for the present 
instrument and evaluated at the various stages of construction and 
modification : 

(i) The ability to make quantitative measurements of both the 
scattered intensity and the temporal intensity autocorrelation 
function for scattering angles ranging from a few mrad down to 
at least 50 /xrad. (We will, in general, specify angular deflections 
in ^rad ; Table I lists conversion factors to other common units 
of angular measure.) 

(it) A stray-light level per coherence solid angle in the scattered 
field (dV.t/dticoH) that was less than 10~ 6 of the incident beam 
power. 
(in) An angular instrumental resolution of less than 15 /urad. 

(iv) The capability of continuously scanning the instrument over 
a reasonable range in scattering angle without the need for 
realignment, 
(y) The attainment of near-diffraction-limited performance using 
customary spherical optics of reasonable cost. 

Taken individually, each of the above goals can be met or bettered 
by existing optical instruments. To cite just two examples, the 200-inch 
Mount Palomar telescope has a diffraction-limit angular resolution of 
about 0.1 /xrad; and, in a typical 6 = 90 degrees light-scattering ex- 
periment, the desired stray-light level would be considered a straight- 
forward achievement. Insofar as the angular range is concerned, we 
can easily show that the scattered light observed at these angles is 
contributed by plane-wave components of the refractive-index per- 



Table I — Conversion factors between various common units 
of angular measure 



1 Deg 
IRad 
1 mrad 
1 /*rad 
1 arc min. 
1 arc sec. 



Deg 



1 

5.73 

0.0573 

5.73 X 10" 
1/60 
1/3600 



Rad 



0.0174 

1 

10~ 3 

io-« 

2.91 X 10"« 
4.85 X 10-« 



mrad 



17.45 

10' 

1 

io-» 

0.291 

4.85 X 10~ s 



Mrad 



17,453 
10« 

10 3 

1 

291 

4.848 



Arc 
Min. 



60 
3438 
3.438 
0.0034 

1 
1/60 



Arc Sec 



3600 

2.06 X 10 8 

206.3 

0.2063 

60 

1 
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turbations in the sample whose wavelengths, A, are given by the 
small-angle Bragg condition 

A = Xo/0, (1) 

where X is the incident-beam wavelength. Therefore, probing the 
scattering-angle range from 50 /*rad to 3 mrad gives information about 
Fourier components of the refractive index having wavelengths be- 
tween 1.0 cm and 0.016 cm, respectively. Here we can point out that 
this spatial-frequency region is routinely examined by common inter- 
ferometric checking methods and holographic techniques. 

The instrument described in this paper is unique in that it meets 
all of the performance criteria simultaneously. In being able to probe 
perturbations with wavelengths as long as 1.0 cm, it represents a 
100-fold improvement on previous low-stray-light-level scattering in- 
strumentation, while its small stray-light background gives it a 
1000-fold sensitivity advantage over conventional interferometric and 
holographic equipment. On a per-unit-aperture-size basis, its ability 
to resolve closely spaced faint (10~ 6 ) and strong (1) features is about 
50 times better than the Mount Palomar telescope. 

2.3 Optical components and physical configuration ot the instrument 

Figure 1 sketches the optical configuration of the most recent 
version of the apparatus designed to meet the performance criteria set 
out in the preceding paragraphs. For brevity, we refer to this par- 
ticular optical system as the MK VI instrument. 




Fig. 1— Optical component layout of the MK VI small-angle-scattering instru- 
ment. Component sizes and spacings are shown approximately to scale. 
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Before discussing the specific function of the various elements of 
the spectrometer, we present below, for reference purposes, a brief 
description of each of these elements and their mounting following the 
identification scheme used in Fig. 1. Whenever spatial or angular 
displacements are specified, they are to be interpreted according to the 
conventions illustrated in Fig. 2. The £ (or 6) and # (or ip) axes are 
taken to be mutually orthogonal cartesian (angular) coordinates 
perpendicular to the axial ray at the point in question. The £(6) 
direction will always lie in the plane of Fig. 1, the instrument's tan- 
gential plane, while $(<p) will denote the vertical or sagittal plane. 
The direction of beam travel defines the local z axis. The basic hardware 
components of the MK VI instrument are the following : 

Ay — An adjustable circular diaphragm stop with an aperture 
diameter cLa\ ~ 5 mm. 

A 2 — A fixed, precision-pinhole aperture, d^2 = 100 /xm. A 2 is 
mounted with £ and # vernier adjustments relative to L\. 

Argon ion laser — The laser is normally adjusted to provide be- 
tween 50 mW and 200 mW of output at either X = 5145 A or 
X = 5017 A. The laser used has a flat-long radius spherical 
resonator, placed at about \ hemispherical spacing, and oscillates 
in TEM o modes only. The output is a well-collimated beam with 
a slight spheroidal distortion. The beam has a gaussian intensity 
profile with a diameter of 1.4 mm as measured to the 1/e 2 points. 




Fig. 2 — Cartesian coordinate system for (x, y, z) showing the cartesian angular 
deflection 6 and <p. 
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j)ET — A silicon-diode photodetector. The diode used is a photo- 
voltaic device operated without bias as a current source. The 
detector has a 1-cm-diameter active area. 

KE — A precision knife-edge custom-fabricated from neutral-density 
"black glass" plate. The 40-mm-long edge is straight to within 
1 fun and nick free. 

Li — Piano convex achromat with a focal length fn = 132 mm. L\ 
and A 2 share a common mount with vernier $ and i degrees of 
freedom. 

jj % — An anastigmatically mounted pair of piano convex achromats 
with an effective focal length fa = 94.77 mm. L 2 has (x , #, i, 0, 0) 
vernier adjustability. 

L 3 — An achromatic lens having fa = 150 mm and a mounted free 
aperture diameter of 35 mm. 

Mi, M%, Ms, Ma — Flat mirrors l\ inches in diameter with X/10 
surface figure. 

Mb, M fl — Dielectrically coated, concave, spherical mirrors fabricated 
of fused quartz. They have a radius of curvature of 2 m and a 
surface conformity of X/10. The mounted free aperture is 6.5 cm 
in diameter. 

Mi, M 8 — Aluminized, first-surface, fused-quartz, flat mirrors. They 
have a mounted free aperture of 13 cm and a surface figure of X/20. 

S — A bilateral slit with straight jaws that can be used to reduce the 
# dimension of the probe beam. 

SL — A commercial, precision, bilateral slit. The jaws have a 50 mm 
usable height and an accurately adjustable opening range from 
3 fim to 3 mm. The slit assembly is mounted on a precision x-z 
translational stage positioned by large-barrel micrometer heads 
with a maximum conforming error of about 1 /mi. The z(0)-axis 
micrometer can be manually positioned or can be driven by a 
digitally controlled stepping motor. 

SV — The location of the scattering sample. 

These optical components are mounted on a 3-inch-thick aluminum 
slab that forms a stable base for the instrument. Because random 
laboratory air currents and temperature gradients can cause angular 
beam deflections comparable to the instrumental resolution, the entire 
apparatus is covered by an essentially air-tight Plexiglas* enclosure. 

2.4 Functional description ot the apparatus 

We can most easily describe the basic optical characteristics of the 
instrument by following the beam path through the system starting 
at the laser source. 



* Registered trademark of Rohn & Haas Company. 
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Mirrors Mi, M 2 , M 3 , and A/ 4 steer the laser output beam around to 
a spatial filter assembly comprised of A\, Li, and Ai. Lens hi and pin- 
hole A 2 form the conventional spatial filter arrangement, while the 
pre-aperture A i serves to block high-angle-beam trash, such as multiple 
reflections in the laser resonator output mirror. The long path length 
through Mi-Mi and aperture A i also provides a significant reduction 
in laser tube discharge light that would otherwise pass through the 
system. 

The spherically spreading wave coming from A 2 is recollimated off- 
axis by Mf,. The collimated beam leaving M b has a diameter* D(l/e) 
of approximately 1.65 cm. The wave-front planarity of this beam is 
measured and adjusted using a wave-front shearing interferometer 
aligned to give a 7-mm shear in the tangential plane. The tangential 
direction wave-front curvature is reduced to less than X/8 over the 
beam aperture by translating the spatial filter assembly along the 
laser beam (I) axis. It is important to note that the use of this off-axis 
collimation scheme produces a large amount of astigmatism and 
tangential plane coma. As a result, it is not possible to make the probe- 
beam wave fronts straight in both the x and # directions simultaneously.* 
The alignment procedure just described is intended to give diffraction- 
limited angular resolution in the 6 plane with some sacrifice in 
direction resolution. 

The collimated probe beam is now sent to the scattering object at 
SV via the flat mirror Mi. Flat mirror M s collects the transmitted 
beam and small-angle scattered light and directs them to Me- 

In the tangential focal plane of spherical mirror M«, the directly 
transmitted beam is brought to a vertical (y) line focus at a position 
we define as xke = 0. Light that has been scattered by some angle 6 
is brought to line focus in the same plane, but at a displaced transverse 
position 

x K b(6) = fm tan d « f M A (2) 

where fm = 100 cm is the focal length of M 6 . Therefore, for sufficiently 
small values of 0, where tan 8 tt d, angular deflection maps linearly 
into lateral displacement at the focus with a position-angle dispersion 
(PAD) constant given by 

PAD(KE) = ^4^ -/*•-! MnV/xrad. (3) 

a 

The knife-edge KE is located in this tangential focal plane with its 
edge vertical and can be set to intercept the transmitted beam at 



* See Section 3.1 for the definition of these quantities. 
f See Section 3.2. 
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6 = to prevent it from entering the remaining portion of the optical 
system. In normal practice, KE is adjusted to occult all light for which 
6 ^ 30 to 50 Mrad. The position and orientation of KE relative to M$ 
is fixed with diffraction-limited accuracy using the standard Foucault 
knife-edge test procedure. 

Lens pair L 2 re-images the focal plane of M e onto the vertically 
oriented main receiving slit SL with a magnification of about (2.54). 
Therefore, the dispersion constant in the slit plane has the value 

PAD(SL) = 2.54MmAxrad (4) 

or 

PAD(SL) = 0.0001 in.//*rad. (5) 

The magnification by L 2 allows the scattering angle 9 to be read directly 
on the "english-units" micrometer that positions the slit. More im- 
portantly, it relaxes the stability and accuracy requirements that must 
be imposed on the slit scan mechanism. Since SL has a minimum 
opening setting of roughly 3 nm, the slit-limited angular resolution is 
about 1 Mrad. 

The proper locations and orientations for the main slit SL and lens 
Lz are determined by an iterative procedure in which one of the jaws 
of SL and the image of the knife-edge formed by L 2 at the slit plane are 
positioned to form an apparent two-jawed slit. The absence of distor- 
tion in the Fraunhofer diffraction pattern formed when this "slit" is 
illuminated by a collimated beam becomes a diffraction-limited test 
for correct lens and slit alignment. 

The scattered light passed by the main slit is collected by L% and 
sent to the photodiode BET. The focal length and position of L 3 are 
chosen such that the real image of the limiting aperture of M 6 formed 
at the plane Ma by lens L 2 is re-imaged onto the detectors active area. 

2.5 Observed angular-resolution performance 

We can assess the 6 direction angular resolution of the MK VI 
apparatus from measurements of intensity as a function of slit position 
(x) in the absence of a scattering object. Two such "instrumental 
profiles" are shown in Fig. 3. The ordinate scale is logarithmic in the 
detector photocurrent with a rough correspondence of 200 /iA/mW of 
optical power. Curves A and B were taken under identical conditions 
except for the position of the knife-edge KE. For curve A, the knife- 
edge was withdrawn to allow the direct probe beam to reach the 
scanning slit, while for curve B, it was positioned to occult all light 
in the region 6 > 50 jurad. Note that the use of the knife-edge provides 
a significant decrease in observed stray-light level, the reduction 
amounting to about an order of magnitude improvement for 6 > 600 
[ira,d (see Section IV). 

1232 THE BELL SYSTEM TECHNICAL JOURNAL, NOVEMBER 1976 




10" 



10" 



10" 



10" 



10" 



10" 



10" 



KT 10 

io- 11 



10 



-12 



1200 



1000 



800 



600 



400 



-200 



6 IN /xRADIANS 



Fig. 3 — Observed instrumental profiles for the MK VI apparatus plotted on a 
logarithmic intensity scale. Curve B was obtained with the knife-edge, KE, occulting 
the direct beam, while curve A was measured with KE retracted. 

An expanded view of the = region of Fig. 3 is shown in Fig. 4. The 
dashed curve represents a best fit of the gaussian exp — [0 2 /£0 2 (l/e)] 
to the instrumental line shape, as detailed in Section 3.1. The full 
width at half-maximum of the fitted curve is 

A0(f ) = 16 M rad. (6) 

For the traces shown in Figs. 3 and 4, the main slit width was set at 

5 Mm which, from eq. (5), is equivalent to a 2-/irad acceptance angle. 

Under these conditions, the effect of artificial slit broadening on the 

line shape may be neglected, as outlined in Appendix A. 

Deriving a value for the sagittal, or <p direction, resolution is a more 

complicated procedure because of the large instrumental astigmatism 

(see Section 3.2). However, a pragmatic number can be given using 

the following operational definition. If the sagittal resolution were 

measured in the tangential focal plane of M 6 , the location of the main 

slit, the instrumental profile would have a full width at half maximum 

given by 

A<p(%) = 93/zrad. (7) 

(See Section 3.2, especially Figs. 16 and 17.) 

The overall angular resolution characteristics of the instrument are 
illustrated in Fig. 5. This sketch shows various fraction-of-maximum- 
intensity contours for the instrumental profile as determined at the 
main slit plane. 
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Fig. 4 — Measured instrumental profiles for the MK VI apparatus in the region 
around = 0. Curve B was obtained with the knife-edge, KE, occulting the direct 
beam, while curve A was measured with KE retracted. Dotted curve is best fit of 
the function exp[— 2 /<50 2 (l/e)] to the transmitted beam profile. 

The measured profiles presented in Figs. 3 and 4 and the corre- 
sponding contours of Fig. 5 were obtained using the full #-axis beam 
height of the instrument, that is, in the absence of aperturing of the 
probe beam by slit S of Fig. 1. As such, the quoted A<p resolution does 
not include any diffraction broadening associated with y direction 
vignetting of the main beam. At full aperture, the instrument's 6 
resolution is essentially diffraction limited, while the q> resolution is 
dominated by astigmatic blurring. However, as the beam height is 
stopped down, diffraction spreading will eventually override the 
astigmatism and the instrument will be solely diffraction limited. For 
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the MK VI, this crossover point occurs at a beam height of about 
0.5 cm or roughly to of the full design aperture. Therefore, the most 
advantageous use can be made of the present apparatus when the 
desired probe-beam geometry consists of a collimated "sheet" or ribbon 
illumination. 

2.6 Analysis of stray-light performance 

A second crucial performance characteristic of any light-scattering 
instrument is its stray-light level in relation to the scattering efficiency 
of the sample under investigation. For the MK VI instrument, the 
ratio of recorded stray-light photocurrent to the photocurrent observed 
at the peak of the transmitted beam, say, can be read directly from 
Figs. 3b and 4b; however, this ratio is not of immediate physical 

<£-/xRADIANS 

A 



1/10 INTENSITY 



1/e INTENSITY 



1/2 INTENSITY 




/^RADIANS 



Fig. 5 — Contours of constant intensity for the instrumental line shape. Each con- 
tour is labelled in terms of a fraction of the peak intensity, I (8 = 0, <p — 0). 
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significance because of the effects of residual instrumental astigmatism. 
The directly measured i(d)/i(0) ratio requires a certain amount of 
mathematical interpretation to provide the stray-light ratio values 
that will be relevant in signal-to-noise or observability calculations. 1-8,27 
In general, the quantities that are most important in this regard are : 
(i) The scattered power per unit solid angle divided by the incident 
probe-beam power : 

JL^ d(P a (d, <p) ^ d(R s (0, <p) (8) 

(Po dddtp dtt 

The ratio [d(R a (d, <p)/d£f\ is a frequently used measure of the scattering 
power of an object; calculating a value of this ratio appropriate to the 
stray light [d(R, ( (0, <p)/dti] provides a basis for estimating the ob- 
servability of a particular scattering feature. This quantity can be 
extracted more or less directly from i(d)/i(0) given (1) the effective 
solid angle subtended by the main slit and detection optics, and (2) 
the instrumental profile contours of Fig. 5. 

(ii) The normalized scattered power per scattering normal mode : 

*ffi>. (9) 

This latter quantity appears in scattered intensity calculations in 
which the index-of -refraction perturbations in the illuminated scatter- 
ing volume are represented in terms of an orthonormal plane wave 
Fourier expansion. 3,27 The mean-square amplitude of these modes and 
their scattering efficiency are, in general, easily calculated from the 
known physical properties of the sample. Characterizing the stray-light 
via a ratio (P, ( (K/)/(Po provides another convenient way of determining 
the observability of the scattering from a particular sample object. 
(Hi) The normalized scattered power per coherence solid angle in 
the scattered field : 

1 d(P >( d > g) = gife g> sw. (10) 

(Po d&coH d£l 

This quantity appears in signal-to-noise ratio calculations relevant 
to determining the spectrum of the scattered light from measurements 
of the temporal autocorrelation function or spectrum of the detected 
photocurrent. 2,3,27 In this case, it is useful to also characterize the 
stray-light level in terms of the quantity [d(R at (d, <p)/(KI C oh2- The 
ratio [d(R„/d£2coff] differs from that defined in eq. (8) in that the solid 
angle is specified as being the solid angle of spatial coherence in the 
scattered field, Q, C oh. The coherence solid angle is a measure of the 
range in 8 and <p about some arbitrary reference direction (6, <p) over 
which the amplitude and/or phase of the scattered electric field exhibits 
statistically correlated behavior. In the typical light-scattering experi- 
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ment, the extent of spatial coherence in the scattered field is controlled 
by the geometry of the scattering sample and wave diffraction. In a 
first approximation, 27 the coherence solid angle is just the diffraction 
solid angle of the scattering source ; that is, 

Qcoh = Adcoii X A<pcoh= (fc-)( 5-) i ( n ) 

where b g and b v are extremal dimensions of the illuminated sample 
volume as viewed from a direction specified by (6, <p). However, for 
the MK VI instrument, the extent of the spatial coherence is partly 
determined by residual aberration effects, and an evaluation of the 
ratio [d(R,t(6, <p)/dQ,coH~] requires a specific calculation of the spatial- 
coherence properties of the optical field at the main slit plane. 

Obtaining expressions for the various stray-light ratios when aber- 
rations are present requires a rather lengthy detailed analysis, as is 
carried out in Section 3.3. For our purposes here, we merely quote 
those results that lead to the numerical ratios appropriate to the 
MK VI instrument. In each case, the procedure is to treat the observed 
stray-light level as if it originated from a fictitious "sample" placed 
at the normal position of the scattering volume. After deriving the 
expressions that relate slit-plane intensity to a real sample's scattering 
cross section, expressed for example as [d<R„(0, <p)/dSf\, we utilize these 
results in reverse fashion to calculate the effective cross section of our 
fictitious stray-light sample. Of course, these expressions derive from 
the main slit-plane imaging characteristics of the instrument; there- 
fore, in succeeding paragraphs, whenever angles or solid angles are 
specified, they are to be interpreted as slit-plane coordinates or areas 
converted to angular units via eq. (5). Consider first the quantity 

d(R 8 t(d, <p) 1 d(? 8l (d, if) 



d&coH (Po ddd<p 



fir 



~ Vo ddd<p McohLvcoh, (12) 

where ABcoh and A<pcoh are the full- width coherence angles in the 
and <p directions. When the main-slit acceptance angles A6sl and A (psl 
satisfy the inequalities 

Ad S L<£ ABcqh = 24.1 /xrad .^ 

A<p S L » A<pcoh = 140.4 /xrad, 

as they do for the profiles of interest here, the right-hand side of eq. 
(12) can be expressed in terms of the measured photocurrent, i{6), as 

d(R 8t (d, ^ = 0) _ iifi) - A^cqh , U n 

dSlcon i(0) x A M1 ' K 

where i(0) is the photocurrent observed at the peak of the direct 
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transmitted beam. In writing eq. (14), we have assumed that A<psl 
is symmetrically placed around <p = 0. The ratio [A <psl/A *pcoh~\ is 
essentially the number of slit-plane coherence areas sampled by the 
detection optics. For the instrumental profiles shown in Figs. (3) and 
(4), A<psl was limited solely by the free aperture of lens L*. Using the 
proper lens free aperture diameter and the linear dispersion constant 
given in eq. (5), we find an effective slit-acceptance angle 

A<Psl = 1.4 X 10 4 /xrad = 0.79°. 

This value of A<psl corresponds to a slit height that samples approxi- 
mately 100 coherence areas. The ratio [_A<psl/ 'A<Pcob1, eq. (14), and 
the data of Fig. 3 combine to give the [d(R st (9, 0)/dQco»] values listed 
in Table II. 

The normalized stray-light power per mode can be found from 
[d(R 8 j(0, 0)/dOcoff] by the methods detailed in Section 3.3. The basic 
procedure involves calculating both the scattered power per coherence 
area and (P 8 (Ky) from a common starting point to obtain the correction 
term that relates them. In the present case, the required relationship 
has the form 



d<R at (d, 0) = (P. t (Kj) A<pcoh ABcoh 
dttcoH <Po (X/&y) (X/&*) ' 



(15) 



where b x and b y are the clear aperture width and height of the instru- 
ment and X is the optical wavelength. The product of the ratios 



AcpcOH 

(X/6 V ) 



and 



Adcc 



(A/6,) 



is a weighted measure of the number of K, modes contributing to the 
power observed in a single coherence area at the main slit. At the full 
aperture of the MK VI instrument, b x = 5 cm and b u = 5 cm, the 
correction factor has the value 



AcpcOH A0COH 

(\/b y ) (X/6.) 



= 33.6. 



(16) 



Table II — Numerical values of various stray-light ratios for the 
MK VI instrument at selected scattering angles 



0-prad 


m/ao) 


da. t (0, 0)/<KIcoh 


<P.,(K,)/<P. 


50 

100 

200 

500 

1000 


1.07 X 10-* 
2.82 X 10~ 5 
9.55 X 10"« 
2.82 X 10" 6 
1.62 X 10~ 6 


1.5 X 10-» 
4.0 X 10~ 7 
1.4 X 10" 7 
4.0 X 10" 8 
2.3 X 10~ 8 


4.5 X 10- 8 
1.2 X 10-i 
4.0 X lO" 9 
1.2 X 10-* 
6.8 X 10" 10 



AdsL = 2/irad, Atp S L = 14,000 ^rad, AOcoh = 24.1 /*rad, and A <pcoh = 140.4 /xrad. 
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Combining eqs. (15) and (16) with the values of [d(R 8 <(0, 0)/dQ C off] 
already calculated gives the stray-light-per-mode ratios to be found 

in Table II. 

The numerical values of the various stray-light ratios may be put 
into perspective by calculating the amplitude of some physical per- 
turbations that would generate a scattered intensity equal to the 
observed stray-light level. Based on a theoretical analysis of the 
scattering problem, it may be shown that, for sufficiently small 
scattering angles, the actual three-dimensional scattering volume can 
be taken to be equivalent to a two-dimensional phase-object placed 
normal to the incoming probe beam. 28 In this two-dimensional phase- 
plate equivalent, the scattering disturbances appear in the form of a 
spatially varying phase thickness ^(x, y), which is the line integral of 
the instantaneous index of refraction encountered by a ray traversing 
the actual sample at the lateral position (x, y). If n(x, y, z) is the local 
index of refraction in the actual three-dimensional scattering problem 
then if/(x, y) is given by 

f(x, y) =~r I ' ™(z, V, z)dz, (17) 

Ao JO 

where L z is the length of the illuminated volume along the direction of 
the incident beam. The phase perturbation fix, y) may be represented 
in terms of a two-dimensional plane-wave Fourier expansion 

Hx, y) = E L U^j)e iK ' x e iK ^, (18) 

Kx Ky 

with the K, = (K x , K v ) chosen to make the expansion functions 
orthonormal over the instrument's full aperture. In this formulation 
of the problem, the normalized scattered power per Ky mode has the 
simple form 

*4^) = i<i^(K,m, as) 

where the angular brackets denote an appropriate time or ensemble 
average. 

The expression for the scattered power given in eq. (19) may be 
used to interpret the stray-light levels observed in the MK VI instru- 
ment in terms of a minimum detectable amplitude for a specific physi- 
cal scattering mechanism. In succeeding paragraphs, we consider three 
such scattering processes : (i) static index of refraction modulation in 
a transparent slab, (it) surface height modulation on a reflecting 
mirror, and (iii) temperature modulation in an otherwise homogeneous 
liquid. 

2.6.1 Refractive modulation in a slab 

The scattering from a static sinusoidal refractive-index modulation 
in a plate is an interesting model problem relevant to holographic 
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memories and general phase-grating problems. We take the sample 
object to be a nominally homogeneous plate of thickness L s and 
refractive index no in which a small sinusoidal index disturbance 

Bn(x, y, z) = Sne iQ " (20) 

has been created. The index perturbation is assumed to be uniform in 
the z direction such that Q lies in the (x, y) plane. The local index in 
the plate is 

n(«,y,i)-n,ri + £**•] (21) 

which, when inserted in eq. (17), gives the phase function \f/(x, y) as 

*(*, V) =Y L [>•>£* + tnLj**} (22) 

A 

The required Fourier amplitude, ( |$(K,) | 2 ), is obtained by inspection 
from eq. (22) as 

fl*™ -[(£)(£ W- (23) 

This result may be used together with eq. (19) to obtain the pertur- 
bation amplitude (8n/n ) necessary to produce a given scattered power 
per mode. For example, taking L z = 1 cm, n = 1.5, and X = 5000 A, 
we find that a refractive-index amplitude dn/no = 1.2 X 10 -9 yields a 
normalized scattered power per mode equal to the MK VI's observed 
stray-light value at 6 = 100 /xrad. To produce scattering at this angle, 
the wavelength of the perturbation A = 2ir/ | Q | would have to be 
A = (Xo/0) = 0.5 cm. Table III lists the "background equivalent" 
bn/ria values corresponding to other values of 0(A). 

2.6.2 Height modulation on a reflecting surface 

Another interesting example from the viewpoint of stray-light level 
comparison is the scattering from a surface height disturbance on an 
otherwise perfect reflecting mirror. Clearly, this problem can also 
serve as a model for calculating the instrumental background when 
mirror surface roughness (see Section IV) is the dominant source of 
stray light. 

Since the primary effect of a surface height deviation is to produce 
a phase perturbation on the reflected wave-front, the phase function 
i[>(x, y) can be written down immediately as 

Hx,y)=^±h{x,y), (24) 

Ao 

where h(x, y) gives the local physical height displacement from the 
nominally perfect geometric surface. 
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Table III — Amplitudes of three scattering perturbations 

necessary to scatter an amount of light equal to 

the instrument's stray-light level 



e 

(/irad) 


A 

(cm) 


<P.(K,)/(Po 


Refractive 

Perturbation 

(Sn/no) 


Surface 

Corrugation 

(dh-k) 


Temperature 
Perturbation 


50 

100 

200 

500 

1000 


1.0 

0.5 

0.25 

0.1 

0.05 


4.5 X 10" 8 
1.2 X 10~ 8 
4.0 X 10-° 
1.2 X 10-» 
6.8 X 10- 10 


2.2 X 10- 9 
1.2 X 10-* 

6.7 X 10- 10 
3.6 X 10" 10 

2.8 X lO" 10 


0.16 

0.084 

0.049 

0.027 

0.021 


6.8 X lO" 8 
3.5 X lO" 4 

2.0 X lO" 6 

1.1 X 10-" 
8.4 X 10" 7 



We will take h(x, y) to be a small, static, sinusoidal corrugation 



h(x, y) = 5he iQ ' T 
for which the phase function is just 



f(x, y) = 



iQt 



(25) 



(26) 



(27) 



The Fourier amplitude <|i£(K y ) I 2 ) follows trivially as 

<|0(K y )|»>= [Y 8h ] 2 ' 

Combining eqs. (27) and (19) with the data of Table II gives the 
"background equivalent" surface corrugation amplitudes listed in 
Table III. Again, these are the surface amplitudes necessary to yield 
a normalized scattered power-per-mode equal to the MK VI's stray- 
light level. For example, when A = 2v/\Q\ = 0.5 cm, the "background 
equivalent" corrugation has an amplitude of 8h = 0.084 A or, in the 
usual surface-figure parlance, 

8h ^ X/60,000. 

2.6.3 Temperature modulation in a liquid 

As a final example, we consider an otherwise homogeneous slab of 
liquid of thickness L z on which is impressed a small sinusoidal tem- 
perature disturbance, 

5T(x, y, z) = We*", (28) 

with Q lying in the (x, y) plane. The calculation of the scattering from 
such an object is really just a simple extension of the result obtained 
above for refractive index modulation. The temperature perturbation 
produces an associated index disturbance that is responsible for the 
scattering. If the temperature perturbation in eq. (28) is impressed 
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isobarically, the associated index modulation is simply 
8n(x, y, z) = ( -^ j ST(x, y, z) 

= ($), m *" (29) 

Equations (20) through (23) may now be used to obtain the relevant 
mean-square phase amplitude, namely, 

<i*«?>iHt„(^),H 2 - (30) 

Taking L. = 1 cm and X = 5000 A and using a typical value for 
(dn/dT) P in liquids, (dn/dT) P =-5 X 10- 4 /°C, we find the back- 
ground-equivalent temperature amplitudes listed in Table III. 

2.7 Conclusion 

In this section, we described the basic features of an optical instru- 
ment capable of extending conventional light-scattering measurements 
to an angular range (50 ,nrad to 3 mrad) not previously accessible. In 
addition to a diffraction-limited angular resolution of a few seconds of 
arc, the MK VI instrument exhibits an exceptionally low stray-light 
background making it an effective tool for probing small-amplitude- 
scattering processes. Besides its primary purpose of studying long- 
wavelength (0.01 cm to 1 cm) thermal fluctuations, the present type 
of apparatus should prove quite useful in other areas where long-wave- 
length perturbations must be probed, such as, 

(*) Holographic and optical memory imaging. 
(tt) Surface roughness testing. 
(Hi) Index of refraction profiling. 

In general, the MK VI offers a sensitivity improvement of a factor 
of about 1000 over the instrumentation normally used for such 
measurements. 

While we have given a rather broad overview of the apparatus in 
the present section, we have not attempted to present the fundamental 
considerations on which the design is based. We refer the reader who 
is interested in these questions to the remaining sections of this paper. 

III. THEORETICAL CONSIDERATIONS IN THE DESIGN OF A 
VERY-SMALL-ANGLE LIGHT-SCATTERING APPARATUS 

3.1 Aperture apodlzatlon 

In the simplest analysis, the ultimate angular resolution of any 
optical instrument is limited solely by diffraction. The expression most 
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widely used to estimate the limiting resolution is the so-called Rayleigh 
criterion, 29 

(A0)(Ad)^X„. (31) 

The quantities Ad and Ad can be interpreted in two ways : 

(i) If Ad is the diameter of a "collimated" beam, then Ad is the 
actual angular spread of the beam imposed by diffraction. 

(ii) If Ad is the collection angle for light emanating from an object, 
then Ad is the smallest spatial detail that can be resolved on 
that object. 

Neglecting for the moment the off-axis features of the actual MK VI 
instrument, we can duplicate its basic function with the two-lens 
system sketched in Fig. 6. Applying the Rayleigh criterion to this 
particular optical arrangement for a collimated beam diameter 
b = Ad = 5 cm, and with X = 5000 A, predicts an instrumental 
angular spread 



A0rayleigh = 10 //rad. 



(32) 



Unfortunately, taken by itself, this value for A0 RAY leigh contributes 
little in the way of a quantitative understanding of the instrument's 
small-angle performance. In fact, the Rayleigh criterion can be mis- 
leading in a number of ways. First, it does not indicate how much light 
an object would have to scatter to be "visible" when the scattering 
angle approaches A0 RAYL eigh- Second, it implies that the angular 
diffraction spread can be decreased to an arbitrarily small value by 
simply increasing the beam diameter b = Ad. In reality, the presence 
of unavoidable optical aberrations will always limit the attainable 
angular resolution. In designing an instrument which is to attain a 
resolution approaching the diffraction limit, a quantitative approach 
to the problem is mandatory. 




Fig. 6 — A simplified "lens equivalent" version of the MK VI optical system. 
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An exact expression for the diffraction-limit resolution charac- 
teristics of an optical system may be obtained as follows. Let us assume 
we know the spatial dependence of the incident beam intensity on 
some surface in the optical system, say the aperture plane of Fig. 6. 
In our case, the electric field on this surface is of the form 

E (x, y) e ^ k ^-^ t \ (33) 

This aperture plane field can be decomposed into a set of infinitely 
extended plane waves, 

E D (d, <p) = E° D (6, p)e-«**"«>«>, (34) 

propagating toward L 2 at various angles with respect to the i axis.* 
The plane-wave amplitudes, E° D (d, tp), are found from the Fourier 
integral 

E° D (0, V) - r f f dxdyEo(x, y)e«°«*+«"«'> (35) 

Ao J JS 

together with the relations, 

| Q | = k = 2ir/X 
Q x = Q . & = k sin 6 at k<0 (36) 

Q y = Q - $ = k smip = k Q <p. 

Equation (35) is just a slightly modified form of the usual scalar 
diffraction theory result which utilizes spherically spreading waves as 
basis functions. 30 

Assume for the moment that lens L 2 in Fig. 6 is infinitely large and 
free of aberrations. Then each of the plane waves, E D (6, <p), is brought 
to a point focus in the slit plane at a position 

£= /tan 0^/0 (37) 

7) = f tan tp ^ ftp, 

where / is the focal length of L 2 , and the approximate signs hold for 
small angles. Combining (35), (36), and (37) gives the field in the 
slit plane as 

Eo(t, n) = 4- [ f E (x, y)e*M*>i*+»>dxdy. (38) 

JAo J JS 

As eq. (38) shows, the field at the slit plane and the field at the aperture 
plane are related as Fourier transform pairs. It should be evident 
that eq. (38) can also be applied "backwards" in Fig. 6 to relate 
E (x, y) to the field at the spatial filter aperture, E (Z', *?')• We will 



* The cartesian coordinate and angle notation follows that adopted in Section II. 
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refer to the field E (x, y) or its intensity 



Io(x, y) = W«o/mo#o(z, y)E* (x, y) 

as the illumination or aperture function. The slit-plane field E (£, y) 
or its intensity, 7o(£, y), is the corresponding instrumental profile. 

The procedure of aperture apodizing may be described simply as 
follows. The basic problem is to find and implement an instrumental 
illumination function, such that both the function itself and its trans- 
form have minimum spatial extent. The goal in a loose sense is to 
optimize the angular resolution per unit aperture opening. Of course, 
one of the general properties of Fourier transform pairs is that the 
second moments or "widths" of the pair members have an approximate 
inverse relationship. The Rayleigh criterion, in fact, is a simplified 
statement of this property. Even within the confines of this inverse 
relationship, however, there is still wide latitude for aperture apodizing, 
i.e., shaping the instrumental profile to obtain particularly desirable 
angular or spatial characteristics. Although the Fourier transform 
relationship between the illumination function and the instrumental 
profile in coherently illuminated optical systems is well known, 31-33 
aperture apodizing schemes are not often applied in optical instrument 
design. Apodizing schemes are, however, extensively employed in 
high-frequency and microwave antenna design, 34,36 where they are 
used to create antenna systems exhibiting an angular directivity 
pattern that satisfies a particular objective. 

In designing an apparatus for very-small-angle light scattering, the 
principal objective is the ability to observe the weak scattered light 
in close angular proximity to the unscattered beam. The goal, then, is 
an instrumental profile that not only has small angular half-power 
points but, more importantly, continues down rapidly to the 10 -5 to 
10 -6 level. The proper shaping of the illumination function, E (x, y) is 
absolutely crucial in obtaining this desired "steep-skirt" behavior. 

In treating the question of aperture apodization for the MK VI 
instrument, we consider the optical system in the simplified form 
shown in Fig. 6. Therefore, the calculated instrumental profiles that 
are obtained below represent the instrument's ideal, diffraction-limited 
performance in the absence of all aberrations. The ways in which the 
residual aberrations of the actual off-axis configuration modify these 
results are taken up in detail in Section 3.2. 

Given the idealized geometry of Fig. 6, the process of evaluating 
various illumination function/instrumental profile combinations can 
be further simplified by the following considerations. First, most of 
the interesting illumination functions and, therefore, their Fourier 
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transforms can be factored to the form 

E (x, y) = E (x)E (y) m) 

Eo(t, n) = Eo(QEo(n). 
Second, in the MK VI apparatus, the open height of the main scanning 
slit guarantees that the measured slit-plane profile is the integral over 
all 7] of the slit-plane intensity Z (£, v)- Under either of these circum- 
stances, we need consider only a one-dimensional form of eq. (38), 
namely, 

Bt(Q = j±p JBt{z)W*+d*, (40) 

where £ = JO. The corresponding aperture and slit-plane intensities are 

2 M Mo ^ (41) 

htt) = I J** ttI [ EoixWO'iWdx 2 - 

^ \ Mo jao| y 

In presenting the results of calculations based on eq. (41), it is 
convenient to adopt a concise terminology to describe the spatial and 
angular widths of the functions involved. We use the following 
notation : 

2>«(£)— Full width at half-maximum intensity for I (x) 
d x {%) — Half width at half-maximum intensity for h(x) 

D f (|) — Full width at half-maximum intensity for 7 (£) or 7 (£') 
d t (i) — Half width at half-maximum intensity for /<,(£) or I (g) 

A0(£)— Full width at half-maximum intensity for Jo(i-) or I (£), 
expressed as an angular equivalent via eq. (37), A0(^) 

= (1//)D«(4) 
86(%)— Half width at half-maximum intensity for / (£) or I (£), 
expressed as an angular equivalent. 

For arguments other than (£) these quantities, give the width at the 
specified fraction of the peak intensity. Since the absolute normaliza- 
tion of the various intensity functions depends only on the total beam 
power, we present all results in terms of the ratio quantities : 

7o(z)//o(0), Jo({)//o(0), /o(0)//o(O), 

where I (6) describes the slit-plane intensity with position £ given in 
terms of the equivalent angular deflection 6 = £//. 

Figure 7 shows the instrumental profiles calculated for two interest- 
ing illumination functions. The first is the uniform field 

^ oW |0; otherwise 
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I (X)/I (0> 
-1.0 



L-j;j_J I L-4. ,1 

-2024 

X-IN cm 



fi IN /xRADIANS 

Fig. 7 — Calculated instrumental profiles for uniform and gaussian illumination 
with b = 5.0 cm and «r chosen to give equal values of &0(\) for both profiles. The 
normalized illumination functions for each case are shown inset on a linear intensity 
scale. 

for which one can easily calculate the following intensity ratios : 

1; -b/2 ^ x ^ 6/2 
; otherwise 
sin 2 (k b!-/2f) 



I (x)/I (0) = 
/o(S)//o(0) = 
/o(0)//o(O) = 



(kM/2f) 2 

sin 2 (k bd/2) 



(42) 



(k bO/2)* ' 

where ko = (27r/X ). The second profile results from a gaussian aperture 
illumination, 

E (x) = Eoe-* 2 ' 2 * 2 

for which the relevant intensity ratios are 



7o(z)// (0) = exp (- x 2 /o*) 
/o(*)//o(0) = exp (- o*%?/f*) 
h(d)/Io(0) = exp (- o*h&). 



(43) 
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The numerical parameters used in obtaining the curves plotted in 
Fig. 7 were 

X = 5000 A, b = 5.0 cm, a = 1.496 cm. (44) 

The b value is representative of the maximum clear aperture of the 
MK VI apparatus ; the value of a was arbitrarily chosen to give equal 
m{\) for both instrumental profiles. Also shown inset in Fig. 7 are 
the two aperture-plane intensity ratios, I (x)/I (0). Note that the 
latter are plotted using a linear ordinate scale. 

As is evident from these two calculated instrumental line shapes, 
the use of gaussian apodization is vastly superior to uniform illumi- 
nation in regard to the observability of weak small-angle features even 
though the two / o (0)//o(O) profiles have identical half- widths. For 
example, at the 10~ 6 level we find 

50(1O -6 ) = 19.8 Mrad gaussian /(x) 
a0(lO -6 ) = 3183 /zrad uniform I(x). 

In fact, from a theoretical standpoint, the gaussian is the ideal form 
of aperture functional. Among the families of possible illumination 
functions, it possesses a unique combination of two properties: (i) 
it has an extremely rapid shirt fall-off, and (it) it goes over into itself 
under the Fourier transform operation. In a general situation where the 
available aperture illumination has some arbitrary (x, y) behavior, 
gaussian apodization would have to be accomplished by interposing a 
suitable neutral density mask at the aperture plane. Fortunately, laser 
sources with a reasonable cavity configuration and oscillating only on 
TEMoo modes have an output beam intensity pattern which is accu- 
rately gaussian, except in the extreme tails of the profile. The avail- 
ability of such a source represents a crucial factor in the feasibility of 
constructing an instrument having the resolution and stray-light per- 
formance of the MK VI apparatus. 

In the actual MK VI instrument, the ratio of the focal length of 
lens L\ to that of mirror M B was chosen to generate a gaussian illumi- 
nation function with an effective width 

<r* = d x (l/e) = 0.826 cm (45) 

in the collimated beam portion of the apparatus (see Fig. 1). The 
instrumental profile calculated via eq. (41) for this value of a is plotted 
in Fig. 8. Also shown for comparison purposes is the profile to be 
expected if we uniformly illuminated the instrument's maximum 
design aperture 

6* = 5.0 cm. 

1248 THE BELL SYSTEM TECHNICAL JOURNAL, NOVEMBER 1976 



It is clear that for these specific values of a and b, gaussian apodizing 
no longer exhibits an absolute superiority over uniform illumination. 
Although the "gaussian" profile still reaches the 10 -6 level much more 
rapidly, it does sacrifice resolving power to the "uniform" profile down 
to approximately the 10 -2 level. What this means to vsa scattering 
performance is the following. The "gaussian" instrument will excel in 
its ability to detect small amounts of light at very small scattering 
angles; however, it will not resolve approximately equal intensity 
features with as much detail as would the "uniform" instrument. As 
we see in the following paragraphs, the tradeoff, roughly speaking, 
involves paying for small-angle weak-intensity performance by sacri- 
ficing some ability to resolve the angular dependent features of the 
scattered light. This comparison can be made more quantitative by 
reference to Table IV, which gives various half-width angles for the 
profiles of Fig. 8. 




60 80 

ft IN /xRADIANS 



Fig. 8 — Calculated instrumental profiles for uniform and gaussian illumination 
with b = 5.0 cm and a = 0.826 cm. The b value corresponds to the maximum clear 
aperture of the MK VI instrument, while the a value corresponds to the width of the 
gaussian illumination actually used in the present apparatus. Normalized illumination 
functions for each case are shown inset on a linear intensity scale. 
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Table IV — Calculated Instrumental profile half-widths at various 

fractions of peak Intensity for three different 

illumination functions 





Gaussian 


Uniform 


Truncated Gaussian 




Illumination 


Illumination 


Illumination 


Xo = 5000 A 


d x (l/e) = <r' 


d x (l/e) = 672 


b = b* a = a' 




<r* = 0.826 cm 


ft* = 5.0 cm 


b/a = 6.05 




(/*rad) 


G*rad) 


(prad) 


50(1/2) 


8.0 


4.43 


8.09 


80(l/e) 


9.6 


5.23 


9.69 


68(10-*) 
80(10r*) 


14.6 


7.36 


14.62 


20.7 


31.8 


20.6 


80(10-*) 


25.3 


100.6 


25.5 


so (io-*) 


29.2 


318 


29.7 


80 (10- 6 ) 


32.7 


1006 


32.6 


60(10-*) 


35.8 


3183 


76.8 


80 (IO" 7 ) 


38.7 


10,060 


243 


80 (IO" 8 ) 


41.3 


31,830 


768 



From the inset plots of I(x)/I(0) shown in Fig. 8, it may seem that 
the gaussian illumination profile used in the present apparatus was 
unnecessarily narrowed relative to the instrumental full aperture. This 
is, in fact, not the case. One crucial detail which has been omitted in 
obtaining the results presented in Figs. 7 and 8 is the possible vignetting 
effect of the instrument's maximum aperture. In calculating E (0 for 
the gaussian E (x), for example, the integral in eq. (40) was taken over 
all x, thereby neglecting any apefturing effects that might occur. 

For the actual vsa scattering instrument, which has a fixed maximum 
aperture, b, E (£), and E (x) are related via the finite domain transform 



(/Ao)* J-b/2 



(46) 



Except for a few special cases, an analytical evaluation of this integral 
is not possible, and one must resort to a numerical approach to investi- 
gate various apodizing schemes. For the experimentally relevant case 
of gaussian illumination, eq. (46) becomes 



EM = 



Er 



•6/2 



(/Xo) 



e -x i l2a i +H2rl/\o)(x^ x _ 
-6/2 



(47) 



On completing the square in the exponential and a change of variable, 
we can rewrite this expression in the form 



£ » ( » - us! ^-<— /;; 



l dw. 



(48) 



1250 THE BELL SYSTEM TECHNICAL JOURNAL, NOVEMBER 1976 



where K, w, w+, and W- are denned as follows : 

x . <rK 



K = (2x/X ) j 

_L_ h 



(49) 



2V2<r V2 
b .aK 



The analyticity of exp (-w 2 ) near w = allows the complex plane 
w; integral to be split into two terms, each having the form of an error 
function of complex argument. Tabulated values of this function are 
available in the literature 36 for a restricted range of the parameters 
(b/a) and (aK). 

In searching out an optimum configuration for the MK VI instru- 
ment, it was decidely more convenient to adopt a fully numerical 
approach in evaluating eq. (47). Appendix B outlines the methods 
that were used. The modified instrumental profile calculations were 
carried out for a range of values of the ratio (b/a) with the aperture 
opening, b, held fixed at b = b* = 5.0 cm. 

Figure 9 shows four such profiles plotted in terms of the normalized 
intensity ratio lo(0)/h(0). Also shown are the corresponding aperture 
ratios J (x)//o(0). The curve for (b/a) = 0.01 is essentially equivalent 
to the result obtained above for uniform aperture illumination. The 
most striking feature of the remaining three h(0)/h(0) curves is the 
presence of an effective background or floor contribution to the profile 
caused by edge diffraction at the aperture. This "shelf" or wing on the 
profile has the slow oscillatory decay of a (sin 2 x)/x 2 functional de- 
pendence. In each case, however, the dttO portion of the curves 
closely approximates the gaussian profile expected from unapertured 
gaussian illumination. 

The results given in Fig. 9 clearly illustrate the tradeoff involved 
in selecting a value of a. In circumstances requiring an instrumental 
line shape with a very low background level, we are forced to accept a 
moderate increase in 50(|) and, therefore, a loss in angular resolving 
power. The curve given in Fig. 9 for (b/a) = 6.05 corresponds to the 
choice that was made for the MK VI apparatus. Various half-width 
values for this profile have been included in Table IV for comparison 
with the results for unapertured gaussian and uniform illumination. 
In the actual instrument, this choice for (b/a) guarantees that the 
calculated edge-diffraction "floor" constitutes less than 10 percent of 
the overall stray-light level. This point is illustrated in Fig. 10 which 
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40 60 80 100 120 

IN juRADIANS 

Fig. 9 — Calculated instrumental profiles for truncated gaussian illumination and 
various values of {b/a). Also shown, on a logarithmic intensity scale, are the cor- 
responding normalized illumination functions \Jo(x)/Io(0)~\. The curves labelled 
(b/a) = 6.05 are appropriate to the b and a values used in the present apparatus. 

shows the theoretical profile for (b/a) = 6.05 superimposed on the 

measured profiles of the MK IV instrument. 

At this point, it is crucial to realize that the truncated transform 

results apply not only downstream from the aperture plane of Fig. 6 

but also upstream toward the spatial filter. There are, in fact, two 

other possible sources of beam vignetting in the system. The most 

obvious is the spatial filter itself. Since the field at the spatial filter 

and the field at the aperture plane are Fourier transform pairs, the 

same considerations involved in choosing (b/a) also apply to the choice 

of spatial filter pinhole size. If edge-diffraction effects at the aperture 

plane are to dominate the system profile, then the ratio of pinhole 

diameter, b PH , to the gaussian focal width at the pinhole, df(l/e), 

must exceed (b/a). Specifically, for the present instrument, b PH must 

satisfy the inequality 

bpH .6 fiftK 

__ > - = 6.05 /K _. 

9.6 Mm a (50) 

bpH > 58/xm. 
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The pinhole diameter actually used is bp H = 100 jum. We note that a 
[bpH/df(l/e)2 ratio this large is contrary to usual spatial filtering 
practices. 

A much more subtle source of possible beam aperturing is the internal 
cavity configuration of the laser source itself. Clearly, the ratio of 
laser tube inside diameter, 6laser, to the mode (1/e) radius, <r LA8ER , 
must also satisfy the inequality 



Olaser ^ 

°LASER a 



(51) 



The laser used in the MK VI has a gaussian mode diameter given as 

^LA8ER(l/e 2 ) = 1.4 mm, (52) 

ii- 




-200 



-100 



400 



I) IN /xRADIANS 



Fig. 10 — Calculated instrumental profile for truncated gaussian illumination 
(dashed curve) superimposed on the measured profiles of the MK VI instrument. 
Curves A and B correspond to the two measurements described in Fig. 4. 
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which is equivalent to 

t«. = Dlu %£ / * ) = 0.5 mm. (53) 

The plasma tube id for the laser is &las E r = 5.0 mm yielding a b/<r 
ratio 

&LASER ^ 1Q (54) 

"laser 

Therefore, in the MK VI instrument, the maximum design aperture 
available to the collimated beam is, in fact, the principal source of 
truncation effects. 

In situations where it is advantageous to alter this resolution-back- 
ground tradeoff by varying (b/<r), it soon becomes apparent that the 
numerical profile calculations of Appendix B are a rather unwieldy 
design tool. Instead, based on an examination of the results shown in 
Fig. 9, it seemed tempting to fit the profile tails to the form 

I T (0) . sin 2 (fcobfl/2) ,_-. 

7W~ l ikjbo/2)* (55; 

and look for an interpolation formula relating the amplitude Ai to 
the ratio (b/a). By a trial-and-error procedure, the following relation 
was found to reproduce the best-fit A i values to within 10-percent error : 



a; 



= i 1 + ST 2 ) ^ ' (56) 

Table V gives the fitted and interpolated values of Ax corresponding 
to the four (b/o) ratios of Fig. 9. It is interesting to note that the 
exponential factor exp ( — 6 2 /4o- 2 ), which dominates the (b/a) de- 
pendence, is just the normalized aperture illumination at the aper- 
ture edge. 

3.2 Optical aberrations 

The fundamental diffraction limitations set out in Section 3.1 are 
really only a prediction regarding the ideal performance of an optical 
system. In the final analysis, the inherent optical aberrations of any 
particular apparatus design determine how close one will come to 
achieving the ideal of diffraction-limited performance. In this section, 
we give a brief summary of those aspects of optical aberration theory 37 
that are relevant to the design of the MK VI apparatus. From a 
qualitative understanding of and analytical expressions for each of the 
various aberrations, we then determine the extent to which aberrations 
modify the ideal diffraction-limited characteristics of the instrument. 
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Table V — Fractional amplitude of the edge diffraction contribu- 
tion to the instrumental profile of truncated gaussian illumi- 
nation. The first column gives the A x values obtained by 
fitting eq. (55) to the tails of the profiles shown in 
Fig. 9; the second column gives the A x value 
predicted by the interpolation formula 
in eq. (56) 



(Wa) 


Ai-Fitted 


A x-Interpolated 



3.333 
6.05 
8.333 


1.0 

1.33 X 10" 1 
5.823 X 10-* 
2.782 X lO" 7 


1.0 

1.49 X 10" 1 
5.92 X 10"* 
2.79 X 10" 7 



Finally, a number of measurements taken on the MK VI instrument 
are compared to the quantitative predictions of the aberration theory. 
Because the present instrument is illuminated with monochromatic 
light, the various chromatic aberrations are absent, and the lowest- 
order non-zero distortions come from the third-order or primary 
aberrations. Here we follow the order-naming convention associated 
with the Taylor expansion of the function sin ^, i.e., 



^3 J,6 



(57) 



where \f/ is the angle of incidence of a ray on a reflecting or refracting 
surface. The approximation sin \f/ = \f/ leads to the usual paraxial optics 
formulae. The next term in the expansion, proportional to ^ 3 , describes 
the primary aberrations. 

The principle aberration-producing elements of the MK VI ap- 
paratus are the off-axis spherical mirrors M 6 and M 6 (see Fig. 1). 
Figure 11 shows the basic optical configuration in which the mirrors 
are used. The labelled geometrical parameters are : 



= ^ The half-field angle or off-axis angle 

a — The semi-aperture 

R — The mirror radius of curvature. 



(58) 



The primary aberrations for an off-axis spherical mirror depend 
parametrically on two angles: u p , the half-field angle and (Q/R), the 
semi-aperture angle. The aberrations associated with the various 
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Fig. 11 — Off-axis mirror configuration used in the MK VI instrument showing the 
important geometrical parameters. 

third-order products of these angles are: 

(&/R) 3 — Spherical aberration 
(a/R) 2 u p — Coma ,^s 

(6,/R)ul — Astigmatism 
u% — Distortion. 

3.2.7 Spherical aberration 

Spherical aberration is a longitudinal focussing defect that is present 
even when the off-axis angle goes to zero. Figure 12 sketches the basic 
ray geometry for a spherical mirror exhibiting pure spherical aber- 
ration. As illustrated in the enlarged detail of the sketch, the marginal 
rays of an incoming parallel bundle are brought to a focus at a point 
closer to the mirror's surface than those lying nearer the axial ray. 
This constantly changing longitudinal focal position results in a trans- 
versely smeared focal spot rather than a focal point. 

One common measure of the amount of spherical aberration is the 
minimum beam waist size produced in the focal region. For a spherical 
reflector, the diameter of this blur spot is given by 

2TSC* = a 3 /R 2 . (60) 

Since transverse displacement at the focus is equivalent to an angular 
deviation in the parallel bundle, we can also express the spherical 
aberration in terms of a full-width angular blur, 

AflSc = —f^- = 2(a/RY, (61) 

where / = R/2 is the focal length of the reflector. 

3.2.2 Astigmatism 

Astigmatism, like spherical aberration, is the result of a longitudinal 
focussing defect. In contrast to the spherical aberration defect, how- 
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ever, the longitudinal focussing error depends not on aperture diameter, 
2($, but on the off-axis angle, u p . Figure 13 illustrates the ray geometry 
of pure astigmatism for a spherical reflector. 

One of the fundamental characteristics of the u p dependent aber- 
rations is the loss of rotational symmetry in the focal region. A non- 
zero off-axis angle destroys this symmetry and establishes two unique 
directions or planes of transverse blurring. The plane defined by the 
incident and reflected axial ray is the tangential plane. Cartesian or 
angular displacements perpendicular to the axial ray and lying in this 
plane are referred to as tangential displacements. The two planes 
orthogonal to this tangential surface and containing either the incident 
or reflected axial ray are called the sagittal planes. Cartesian or angular 
displacements from the axial ray in these planes are sagittal 
displacements. 

For a spherical reflector exhibiting pure astigmatism, a fan of 
parallel tangential plane rays are brought to a focus closer to the mirror 
surface than an identical sagittal fan. The focal region pattern found 
by decomposing the entire illuminated aperture into such ray fans 
consists of the two longitudinally separated focal lines depicted in 
Fig. 13a. The longitudinal (2) separation of the *S and T foci (2 AC*) 



PARAXIAL 
RAYS 




/ SPHERICAL 

/ ABERRATION 



PARAXIAL 
FOCUS 



/ MARGINAL 
FOCUS 



Fig. 12 — Ray diagram for a single spherical mirror exhibiting a pure spherical 
aberration defect. 
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can be calculated from the relation 

AC* = ^p (62) 

and depends only on the semi-field angle. 

In the absence of other aberrations (and the effects of diffraction), 
the S and T focal lines are infinitely narrow in their respective planes. 
This means that the sagittal height of the tangential focus and the 
tangential width of the sagittal focus can be obtained from simple 
extremal ray geometry. For example, Fig. 13b shows the extremal rays 
seen in a tangential plane projection. Normally 2AC* is small compared 
to the reflectors focal length (R/2). It follows that the lengths of the 
two focal lines are identical and given by 

2TAC* = 2AC* (|^) = 2ai£ (63) 

The full-width angular spread equivalent to this spatial blur is 

Mtac* = 4(a/R)u*. (64) 

3.2.3 Coma 

When the off-axis angle is non-zero, the longitudinal focussing error 
that produces spherical aberration also gives rise to an asymmetric 
transverse blurring called coma. Figure 14a sketches the basic elements 
of the focal region pattern for a spherical reflector exhibiting a pure 
coma defect. Rays in the paraxial region are brought to a focus at the 
axial focus, P, while rays from larger-diameter annular zones on the 
mirror's surface form focal circles whose centers are tangentially dis- 
placed from P. The radius of a particular focal circle increases as the 
square of the radius of the zone producing it. 

Figure 14b gives a qualitative representation of the characteristics 
of the focal pattern as found by dividing up the illuminated aperture 
into these annular zones. Mathematically speaking, the focal circles 
are not sharp unless the radial thickness of the corresponding zones 
vanish; however, the sketch does predict quite nicely the overall 
exterior outline of the coma blur patch. 

The continuum of focal circles nest into a 60° wedge extending out 
from the axial focus forming a pattern commonly called comatic flare. 
The largest-diameter focal circle, produced by the annular zone at the 
edge of the illuminated aperture, has a radius CC* given by 

CC* = ^2- (65) 



2R 
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(b) 



2 SCC 




Fig. 14 — Partial ray diagram, (A), of an off-axis spherical mirror exhibiting a 
pure coma defect. The details of the focal plane pattern are illustrated in (B). 



Its center is tangentially displaced from the axial focus by an amount 
2CC*. 

Although coma is a highly asymmetric aberration, it is still con- 
venient to specify its effect in terms of transverse and/or angular blur 
sizes. The numbers ordinarily quoted for coma correspond to the 
extremal dimensions of the coma patch in the tangential and sagittal 
directions. It follows easily from Fig. 14b that the full-width trans- 
verse spatial blurs are 



f direction 2TCC* = 3CC* = 
S direction 2SCC* = 2CC* = 



3a 2 u, 
2R 

a 2 u p 
R 



(66) 



The equivalent full-width blur angles are 

T direction Ad* TC c = S(d/R) 2 u p 
S direction A0s CC = 2(a/R) 2 u p . 



(67) 
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3.2.4 Distortion 

The last of the primary aberrations is a defect of off-axis magnifi- 
cation called distortion. This particular aberration is associated princi- 
pally with optical systems that form a real image of an extended object 
at finite magnification. For example, the image of a rectangular grid 
of regularly spaced points will exhibit the classic "barrel" or "pin- 
cushion" appearance in an optical system involving pure distortion. 
For the situation of interest here, namely a spherical reflector with 
the image at infinity or at the focus, the amplitude of the pure distor- 
tion aberration vanishes identically. 

3.2.5 Application of the aberration results to the MK IV instrument 

The third-order aberration theory results for the off-axis spherical 
reflector are summarized in Table VI which gives the expressions for 
the various transverse and angular blurs. It must be emphasized that 
the aberration theory results outlined above are derived from purely 
geometric ray tracing. In no sense does this theory predict the actual 
intensity distribution in the image plane for a specific aperture illumi- 
nation. The transverse and angular blur patterns define the outlines 
of a boundary between focal illumination and strict geometric shadow 
in the absence of all wave interference and diffraction effects. However, 
in certain situations, it is possible to combine the geometric aberration 
results with aberration-free diffraction calculations to obtain useful 
instrumental profile information. The approach works well when one 
or more of the following conditions are satisfied : 

(i) Diffraction blurring is large or small compared to spherical 

aberration. 
(u) One primary aberration is dominant, 
(m) The ideal diffraction-limited system profile is free of large 

interference maxima and minima. 

Table VI — Summary of the analytic expressions for the 

transverse and angular aberration blurs for a single 

off-axis spherical mirror 



Aberration 


Full- Width Transverse Blur 


Full-Width Angular Blur 


Spherical 


a*/R* (2 TSC*) 


2(a/R) 3 


Coma 


t -3a*u P /2R (3CC*) 
S -2Qhi p /2R (2CC) 


-3(a/R)hi p 
-2(a/R)hi P 


Astigmatism 


2du p (2 TACT) 


4(a/K)?4 


Distortion 


BvlR 


2Bvi 
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In these cases, one can intuitively construct the aberration-affected 
profiles with a fair degree of accuracy using the one-to-one geometric 
mapping of regions of the aperture onto the focal plane. For example, 
it can be argued from Fig. 13b that for pure astigmatism the tangential 
direction intensity profile at the sagittal focus should be a demagnified 
replica of the aperture illumination. This type of analysis is essential 
in obtaining quantitative results from the aberration expressions. 

In the MK VI instrument, the off-axis angle and radii of mirrors 
M 6 and M* are 

7 = 2\u p \ = 0.116 rad ,„ Q . 

R = 200 cm. (68) 

Assigning a value to be used for the semi-aperture (J is a more subtle 
question, especially since we are interested in gaussian rather than 
uniform aperture illumination. However, in the spirit of the astigma- 
tism example given in the preceding paragraph, we take 

a = d x (l/e) = a* = 0.826 cm (69) 

and assume that the focal plane profiles will also be gaussian. On the 
basis of geometrical imaging, the blur values calculated via Table VI 
should then be (with the exception of coma) the full-width to the 
(1/e) points of a focal plane gaussian profile. We show in succeeding 
paragraphs that these assumptions lead to a selfconsistent picture 
of the experimentally observed aberration effects in the MK VI 
instrument. 

Table VII gives the transverse and angular blurs for a single spherical 
reflector in the MK VI configuration. The table also includes the 
longitudinal separation of the S and T focal planes as well as the full 
angular width of the focus imposed by diffraction. From the viewpoint 
of small-angle-scattering performance, the most serious of the aber- 

Table VII — Numerical values of the aberration blurs for a single 
spherical mirror used in the MK VI configuration 



Aberration 


Full-Width Transverse 
Blur (/xm) 


Full-Width Angular Blur 
(/*rad) 


Spherical 


0.141 


0.141 


Coma 


f 2.96 
S 1.97 


2.96 
1.97 


Astigmatism 


55.6 


55.6 


Distortion 









OS-rWration = 2 AC* = 0.336 cm. A0(l/e) d iHraction = 19.3 /«rad. 
1262 THE BELL SYSTEM TECHNICAL JOURNAL, NOVEMBER 1976 



Table VIM — Numerical values of the aberration blurs for the 

system of two off-axis spherical mirrors used 

in the MK VI apparatus 



Aberration 


Full-Width Transverse 
Blur (jtva) 


Full-Width Angular Blur 
(/irad) 


Spherical 
Coma 

Astigmatism 
Distortion 


0.282 

T 

S 

111.2 



0.282 




111.2 





OS-TWmtion = 2(2 AC) = 0.672 cm. A0(l/e)dirfr ac uon = 19.3 /irad. 
rations is coma. Even though the calculated comatic blurs are numeri- 
cally small compared to the diffraction spread, the presence of coma 
can result in a distinctly asymmetric instrumental profile. Moreover, 
in a coherently illuminated system, the coma flare is criss-crossed by 
interference patterns whose tails extend far beyond the calculated 
geometric limits. This latter effect can significantly raise the effective 
"floor" level of the instrumental profile. 

Fortunately, in a symmetric two-mirror system like the MK VI, 
the geometry may be chosen such that the total coma vanishes identi- 
cally. In fact, all aberrations that depend on an odd power of the half- 
field angle u p disappear if the field angles at the two elements are 
made equal and of opposite sign. By convention, u p is defined as the 
angle through which the incoming axial ray must be rotated to bring 
it into coincidence with the local radius vector of the element's spheri- 
cal surface (see Fig. 8). An inspection of Fig. 1 shows that in the MK VI 
apparatus the field angle rotations at M 6 and M& are of opposite sense. 
In this case, when the two off-axis angles are made equal in magnitude 
the coma and distortion aberrations vanish while the spherical aber- 
ration and astigmatism double. 

The total calculated aberration blurs for the instrument are sum- 
marized in Table VIII. Since the spherical aberration is small com- 
pared to the diffraction spread, it is reasonable to expect that the 
interpretation of these blur values as (1/e) full- widths of a gaussian 
blur profile should work quite well. This is in fact the case. 

The tangential and sagittal foci of the MK VI instrument were 
located using a modified Foucault knife-edge procedure, and T direc- 
tion scans of the intensity profiles were taken in each case. The mea- 
sured S-T separation was 

(S-T) separation = 0.660 ± 0.013 cm. (70) 
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Fig. 15 — Observed instrumental profile as measured by a tangential direction 
scan in the instrument's tangential focal plane. The heavy dots are a best fit to the 
function exp [-«*/»■ (l/e)j\ 

At the T focus, only diffraction and spherical aberration contribute 
to the profile width. From the standard gaussian convolution formula, 
we can calculate the expected A0(l/e) : 



T PLANE 



A0(l/e) = V(19.30) 2 + (0.282) 2 M rad 



T SCAN 

= 19.302 M rad. (71) 

Clearly the spherical aberration has a negligible effect on the ideal 
diffraction-limited broadening. Figure 15 shows a typical high-resolu- 
tion jT-plane scan for the instrument. The large apparent noise in this 
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trace is produced by residual air currents and vibration in the ap- 
paratus and corresponds to a peak-to-peak beam wander of roughly 
0.5 /xrad (0.7 arc second). The results of fitting the observed profile 
with a gaussian shape are indicated by the points in Fig. 15 and give 
an experimental full width 

A0(l/e) = 250(1 /e) = 19.2 M rad. (72) 

In the S focal plane, diffraction, astigmatism, and spherical aberration 
all contribute to the instrumental line shape. The full-width calculated 
from Table VIII is 



S PLANE 
T SCAN 



A0(l/e) = V(19.3) 2 + (0.282) 2 + (111.2) 2 M rad 

= 113 Mrad. (73) 

The predicted width comes predominately from the astigmatic blur- 
ring. Figures 16 and 17 give two experimental S-plane profiles recorded 
with logarithmic and linear intensity scales, respectively. A gaussian 
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Fig. 16 — Observed instrumental profile as measured by a tangential direction scan 
in the instrument's sagittal focal plane. The heavy dots and dashed curve a are 
best fit to the function exp [— </?/ 6 {/?(!/ eft. 
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Fig. 17 — Observed instrumental profile as measured by a tangential direction scan 
in the instrument's sagittal focal plane. The heavy dots are a best fit to the function 
8xp[-^/M(l/«)]. 

fit to the logarithmic curve, indicated by the points in Fig. 16, gives 
very good agreement with the observed line shape over roughly four 
orders of magnitude in intensity. The range and precision of the fit 
provide strong support to our assumptions regarding the interpretation 
of the aberration blur values. The best-fit half-width values for the 
logarithmic and linear scans are 86 (1 /e) = 57 fir ad and 69(l/e) = 55 
Mrad, respectively. The mean observed full-width 

A6(l/e) = 112 Mrad 

is in excellent agreement with the calculated value given in eq. (73). 
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3.2.6 Conclusions 

In summarizing the discussion of aberrations, a number of general 
points deserve to be made and reiterated regarding the relationship 
between aberrations and small-angle-scattering performance. 

(i) The cancellation of the asymmetric aberrations, coma especially, 
is crucial in obtaining an instrumental profile that has symmetry, the 
necessary steep skirt fall-off, and a low background value. 

(ii) Spherical aberration, although it has a negligible effect relative 
to diffraction in the present instrument, can rapidly grow to serious 
proportions with increasing aperture size, ft = d x (l/e). The angular 
blur of this aberration, A0^ c , increases as the cube of the aperture size 
while the diffraction spread varies inversely with d x {\/e). The relative 
contribution of spherical aberration to the profile width will, therefore, 
increase as (2 4 . Since A0s C depends only on the reduced quantity ( 0,/R), 
however, a constant ratio of spherical blur to diffraction spread can 
always be obtained by scaling the mirror radius R to keep &*/R 3 
constant. For example, an instrument with 10 times better angular 
resolution than the MK VI might conceivably utilize 80-cm-diameter 
mirrors with a 21.5-meter focal length. 

(Hi) The presence of a large residual astigmatism need not be 
detrimental if one is satisfied with an instrumental performance that 
is diffraction limited in only a single angular direction. It might appear 
from Fig. 13 that simultaneous sagittal and tangential resolution could 
be achieved by placing separate slits at the S and T foci of the collecting 
mirror. This is true if the wavefronts of the incoming ray bundle are 
perfectly parallel. However, in a two-mirror symmetric apparatus, such 
a bundle cannot be produced because of the collimating mirror astig- 
matism. For example, with reference to Fig. 1, the spatial filter pinhole, 
A 2 , can be placed at either the T or S focus of mirror M B . In the first 
case, the wavefronts of the beam travelling toward M 1 are tangentially 
collimated but sagittally curved ; in the second case, the converse is 
true. From the viewpoint of light-scattering kinematics, this "col- 
limated" beam will be able to conserve momentum with a relatively 
broad range of scattering vectors lying in the plane containing the 
wavefront curvature. Thus, even though mirror il/e forms S and T 
focal lines of equal sharpness, only a single high-resolution axis actually 
exists for either position of pinhole A 2 . 

In the MK VI instrument, the spatial filter pinhole is at the T 
focus of the collimating mirror so that the probe beam is tangentially 
collimated. The nature of the wavefront curvature in the sagittal plane 
may be calculated in a simple fashion from the known (S-T) separation 
2 AC*. Since the pinhole (at the T focus) is closer to the mirror's surface 
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Fig. 18 — Illustration of the wavefront curvature existing in the off-axis mirror 
collimating arm of the MK VI instrument. 

than the S focus, the sagittal plane wavefronts appear to diverge from 
a point source lying behind the mirror. The object distance, q, between 
this virtual point source and the mirror, may be calculated via the 
usual paraxial formula. Since (S-T) = 2 AC* is small compared to the 
focal length of M B , we have 



Q = 



P 



R* 



2AC* 4(24C*) 



(74) 



Inserting the appropriate numerical values, R = 200 cm and 2 AC* 
= 0.336 cm, gives 



q *£ 3.98 X 10 4 cm. 



(75) 



Figure 18 shows a sketch of the probe-beam constant-phase surfaces 
with the sagittal curvature greatly exaggerated for clarity. The magni- 
tude of the curvature can be specified in terms of the longitudinal 
spatial separation between the wavefront and a reference plane which 
is tangent to the wavefront at the axial ray. Since the virtual sagittal 
source point lies so far behind the mirror, it does not matter exactly 
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where in the optical path we calculate this separation. The maximum 

deviation between the wavefront surface and the reference plane, Vz, 

occurs at the sagittal extremes of the beam and is easily found to be 

given by 

ft 2 
Vz = _. (76) 



If we take d = b*/2, where b* = 5.0 cm is the maximum clear aperture 
diameter in the MK VI instrument, we have 

Vz S* 10.5 X 10~ 5 cm, 

or roughly 2X of peak deviation from perfect collimation. In a situation 
where the full sagittal (y) aperture height need not be used, it is 
possible to reduce this deviation substantially since Vz is propor- 
tional to ft 2 . For example, if we aperture the height of the probe beam 
to b = 5 mm, the peak wave-front deviation is reduced to Vz = X/50, 
or essentially perfect collimation. Thus, the most desirable probe- 
beam configuration in the MK VI instrument corresponds to a "flat 
ribbon" or "sheet" type of illumination. 

3.3 Scattered field intensity, spatial coherence, and scattering kinematics 
In the presence of aberrations 

In a light-scattering optical system whose angular resolution capa- 
bilities are in some respect dominated by aberration effects — for 
example, the astigmatic <p blurring in the present instance — we find 
that other important properties of the observed scattered field are 
modified by the aberrations as well. In this section, we examine three 
aspects of normal light-scattering theory that are qualitatively altered 
by the presence of aberrations : 

(i) The form of the spatial coherence function for the scattered 

field, 
(it) The application of the normal kinematic restrictions (or wave 

vector conservation conditions) in the scattering process. 
(Hi) The calculation of the amplitude of the perturbations that 
give rise to observed levels of scattered or stray light. 

The effects of aberrations in all three cases have a straightforward 
physical interpretation connected with the fact that light scattered 
into a specific direction (6, <p) is no longer brought to a diffraction- 
limited spot focus at the observation plane. 

In typical calculations of the scattered field, in which an incoming 
plane wave is assumed to impinge on the sample, the far-field angular 
distribution of scattered intensity is shown to be simply the spatial 
Fourier transform of the refractive-index perturbations in the illumi- 
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nated volume. 13 The light observed at a particular angular position 
(0, <p) is contributed by a Fourier component of the refractive index K 
given by 

K = k, - k 0> (77) 

where k is the wave vector of the incident beam and k„ points in the 
direction of observation (6, <p) and has a magnitude |k«| = |k | 
= 27r/X . For a finite illuminated volume, the refractive-index per- 
turbations are most usefully represented in terms of a plane-wave 
Fourier expansion 

Z $Uj exp (tK r r'), (78) 

i 

with the Ky chosen to make the expansion functions orthonormal over 
the scattering volume. 27 The scattering of a collimated incident beam 
by this assembly of plane waves consists of a family of diffracted 
beams that originate from the K/s satisfying the Bragg condition, 
eq. (77). On the surface of a sphere in the far field, these diffracted 
beams form a contiguous but essentially nonoverlapping series of 
diffraction "spots," each associated with a particular K,. In the usual 
situation, where the amplitudes of the individual Ky disturbances are 
statistically independent, these patterns also delineate areas or solid 
angles of statistical field correlation. If the far-field scattered radiation 
is focussed onto the observation plane by an ideal lens or mirror, this 
contiguous angular distribution of "spots" is imaged one-for-one onto 
the focal plane. A ray penetrating the reference sphere at an angular 
position (0, <p) is imaged onto the focal plane at a transverse position 
(£> v), where, in the small angle limit, 

* = fO, v = f<P- (79) 

In this ideal situation, the intensity observed at some (£, 17) is scattered 
essentially by a single K,- plane-wave mode. The measured intensity 
may, in theory, be used to calculate the mean-square-amplitude of 
the mode, or vice-versa. Furthermore, the spatial coherence properties 
of the field at the observation plane are determined uniquely by the 
angular distribution of intensity within one of the diffracted beams. 
Formally speaking, the presence of aberrations in the imaging of 
the reference sphere scattered field produces qualitatively the same 
effects as any other imperfect focussing of the far-field pattern. The 
pattern of diffraction spots will be formed with a degree of spot 
broadening and overlap that depends on the nature and extent of the 
focussing defect. The scattered light reaching a specific (£, 77) point 
at the observation plane is no longer associated with a single K ; 
disturbance, but is an appropriately weighted sum of contributions 
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from a number of modes. As a result, the Bragg condition, eq. (77), 
is not strictly applicable in relating a particular (£, 77) observation 
point to a specific plane-wave disturbance. The aberration or defocus- 
sing effect must be understood in detail before the measured intensity, 
or its time evolution may be used to infer the physical behavior of 
modes responsible for the scattering. Under defocussed conditions the 
spatial correlation function is also modified, though its functional form 
does retain a close resemblance to the intensity pattern associated 
with a single Ky diffraction "spot." In the following paragraphs we 
consider how the residual astigmatic blurring in the MK VI instrument 
affects the three slit-plane field properties enumerated in the opening 
paragraph. 

3.3.1 Kinematic relations 

To understand how the wave vector conservation criterion is to be 
applied at the slit plane of the present apparatus, we need to know 
(i) the slit-plane intensity pattern formed by scattering from a single 
plane-wave disturbance, and (u) the relative positioning of the spots 
from the various allowed Ky. In the MK VI instrument, the intensity 
pattern associated with a single Ky is identical to the diffraction- and 
aberration-affected instrumental profile whose properties were treated 
in detail in Sections 3.1 and 3.2. At the slit plane, therefore, the single 
Ky diffraction spots have the elongated gaussian shape depicted in 
Fig. 5. Expressed in terms of angular coordinates via eq. (79), the 
normalized intensity distribution within a "spot" is simply 

m g) _ cxp I (g - W 1 c::p (_ (y-w)M . (80) 

The reference point (0y, <p } ) specifies the angular position of the spot 
center, which, in the present case, is correctly predicted by the Bragg 
condition, eq. (77). 

Given a correct form for the intensity distribution within a single 
Ky pattern, we must still determine the slit-plane spacing of the spots 
associated with the family of allowed Ky. Clearly, this spacing depends 
on the reciprocal lattice of the orthonormal expansion functions 
exp (zKyr') which, in turn, is fixed by the geometry of the scat- 
tering volume. For sufficiently small scattering angles, the actual 
three-dimensional scattering sample can be taken to be equivalent to 
a two-dimensional phase object placed normal to the incoming probe 
beam. 28 The scattering disturbances in this "phase-sheet" may be rep- 
resented in terms of a two-dimensional plane-wave Fourier expansion 

EI • • -exp (iK x x) exp (iK y y), 

K x Ky 
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with the (K x , K y ) chosen to make the expansion functions orthonormal 
over the instrument's full aperture. For a rectangular aperture with 
full-width x and y dimensions b x and b y , the allowed Ky can be obtained 
from the cyclic boundary condition relations 

K x (m) = — , K v (n) = -y- , (81) 

where m and n are the integers 

m, n = 0, ±1, ±2, •••. (82) 

The scattering angles (0/, p,) for the central ray of the diffracted beam 
produced by a particular K(m, n) are then given by the small-angle 
Bragg conditions, 

-rr , ^ i n 2irm 

K x (m) = k(fi m = -r— 

0x (83) 

v ( N , 27rn 

A„(n) = fco^n = -j— ' 

°y 

It follows that the family of diffracted beams are brought to a focus 
at the slit plane on the vertices of a rectangular mesh whose grid 
spacings are given by 

&p = f0 8 p = /(£-°) (84) 

Vsp = f<psp = f f t- J ■ (85) 

For the MK VI instrument at full aperture (b x = 5.0 cm, b v = 5.0 cm), 
the equivalent angular mesh spacings are 

Osp = <psp = 10 Mrad. (86) 

If we imagine the instrumental profile contours of Fig. 5 arranged on 
such a mesh, there will be little overlap in the 6 direction but consider- 
able overlap along 0. The light received at some (6, <p) point in the 
slit plane will contain contributions from roughly 10 distinct K(m, n), 
each having the same m(K x ) index but differing n(K u ) components. 
Because the MK VI instrument is capable of probing values so close 
to the diffraction limit, corresponding to a very small m index, 
m = 5-300, these multiple contributions can prove a serious problem. 
This point is illustrated in Fig. 19, which shows the (1/e) contour of 
a single K(m, n) intensity pattern centered at = 80 ,urad, <p = 
superimposed on the slit-plane mesh of the (0 m , <?„)• The and <p axes 
of the figure can also be labelled in terms of the wave vector com- 
ponents K x and K y to which the angles are directly proportional as 
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Fig. 19 — Resolution function of the MK VI instrument superimposed on the 
slit-plane mesh points corresponding to the allowed scattering vectors K,-. 

indicated in eq. (83). For the eleven mesh points falling within the 
contour, the magnitude of the associated K(m, n) varies considerably. 
For the point at the center of the contour, we have (m = 8, n = 0) and 

|K| = 2ir[(m/b x ) 2 + (n/6„) 2 ]i = 10.0 cm" 1 , 

while for the points at the <p extremes (m = 8, n = 5), we find 

|K| = 11.9 cm- 1 . 

When the physical properties of the modes are strongly | K | dependent, 
this overlap can lead to a difficult task in the interpretation of the 
measured intensity and/or its time dependence. 

Clearly, the | K | smearing effect becomes less significant as 6 in- 
creases. Less obvious is the fact that the problem of multiple K, 
contributions can be alleviated by stopping down the beam height of 
the instrument, b y . As was pointed out in the conclusion of Section 
3.2, the c}iindrical distortion of the probe-beam wave fronts, which is 
a manifestation of the collimating mirror's astigmatism, can be made 
negligibly small by reducing b y . With b y ^ 0.5 cm, for example, the 
beam incident on the sample can be considered as collimated to within 
the diffraction limit. In this case, the usual kinematic conditions apply 
in relation to the far-field scattered light ; that is, the far-field array of 
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"spots" now form a contiguous and nonoverlapping pattern. Of course, 
each spot is elongated along <p by diffraction-spreading because of the 
imposed asymmetry of the probe-beam dimensions, i.e, b x = 5.0 cm and 
b y ^ 0.5 cm, but the angular spot spacings d S p = (X/6 X ) and 
VSP = (\/b y ) are correspondingly asymmetric. Given this particular 
situation in the far field, we must still consider the effect of the astig- 
matism associated with the light-collecting mirror. As b y is decreased, 
the direction astigmatic blurring at the slit plane is decreased pro- 
portionally,* while <p blurring due to diffraction increases. At some 
point, a crossover occurs beyond which diffraction spread dominates 
the slit-plane imaging. In this limit, the collecting mirror appears 
aberration free and the far-field pattern of K, spots undergoes the 
normal one-to-one, no-overlap mapping onto the slit plane. For the 
MK VI configuration, the crossover occurs at b v « 0.5 cm or roughly 
^ of the full design aperture height. Of course, it should be noted that 
under fully diffraction-limited conditions (b x = 5.0 cm, b v < 0.5 cm), 
the instrument retains its very asymmetric resolution profile. What we 
have done is to make the instrument appear to be in "good focus" by 
introducing a sufficient amount of <p direction diffraction spreading to 
swamp the aberration defocussing. The price paid for this is that the 
instrument becomes incapable of probing scattering disturbances 
having as long a y direction wavelength (that is as small a value of K v ) 
as can be resolved in the £ direction. 

3.3.2 Relation between the silt-plane intensity and the scattering 
cross-section ol the sample 

In the absence of aberrations or other defocussing problems, the 
far-field or observation plane scattered intensity can be related easily 
to the mean-square amplitude of the scattering perturbations using the 
standard integral expression for the scattered field. 1 - 3 - 27 In the presence 
of imaging aberrations, the total scattered power per plane-wave mode 
is unchanged ; however, now the calculation of the slit-plane intensity 
is complicated by the overlap of the various K> diffraction spots. For 
the MK VI instrument, a relation between the scattered power per 
mode and the observed slit-plane intensity may be obtained as follows. 

The slit-plane intensity corresponding to a single K y disturbance 
can be written down formally as 

""-'£..(*. *>I-P [- ( j^} exp [- fe^] (87) 

with slit-plane position specified in angular units. The position of the 
* See Figure 13 and the discussion pertaining to this figure. 
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central ray of the pattern (0 m , <p„) is given by eq. (83). The total 
scattered power in the single Ky pattern is just 

P.(Ky) = J Jddd<p\E nm (d, <p)\ 2 

= T5d(l/e)8<p(l/e)\E nm (e m , <p n )\\ (88) 

The problem now is to relate the observed total slit-plane intensity 
to the peak mode intensities, \E„ m (9 m , <^ n )| 2 . The intensity observed 
at the position (6, <?) is the sum over all possible mode contributions : 

= EJ^(^^Pexp[-^^]exp[--^^]. (89) 

If \E„ m (d m , ip n ) | 2 is independent of (n, m) over the range where the 
gaussian terms are nonvanishing, eq. (89) can be simplified to give 

\E.(fi, ^)| 2 = \E nm (6 m , <p n )\*S m S n , (90) 

where S m and S„ are the factored sums 

S m = E exp - {Id - m(\o/b x )J/W(l/e)) (91) 

tn 

S n = E exp - {[> - n(\ /b y )J/?><p*(l/e)}. (92) 

n 

Consider the »S m sum expressed in the following dimensionless form 
= Eexp (-[« - mo] 2 ), (93) 

m 

where 



50(1 /e) ' 50(l/e) 6, 

Although the indicated summation cannot be carried out explicitly, 
it can be expressed in a more useful form via the identity 

E exp [— (u — ma) 2 ] = — E exp ( — J cos (2mmi/a), (95) 

which is an immediate corollary of Poisson's formula. 38 Applying this 
identity to the S m summation yields 



Q V^fl(l/e) ^ 

^ = -ix7^r? exp 



*gg(l/g) 1 noo f 2nnrd "1 

W^~J C0S [(x7^)J" (96) 
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For the numerical parameters relevant to the MK VI instrument, 
60(1 /e) = 9.6/irad, (X/6«) = 10 vr&d, we need retain only the m = 
term in the right-hand side of eq. (96). With an error less than 10~ 4 , 
we may take 

A similar result follows easily for S n . We then have for the total slit- 
plane intensity 

l*W ^ ~ 1^-^-' ^ n) l (x/6 x ) (X/b y ) (98) 

As this result shows, the various elongated gaussian patterns associated 
with the individual K y modes overlap in the slit plane to produce an 
essentially uniform illumination. 

The unknown intensity factors | E nm (6 m , </>„) | 2 can now be eliminated 
between eqs. (88) and (98) to give the desired relationship between 
the overall slit-plane intensity and the scattered power per mode, 
namely, 

P.(Ki) = \E t (d h wWCVWCXA). (") 

Equation (99) is the basic result which allows the measured intensity 
to be related quantitatively to the amplitudes of the individual 
scattering perturbations.* 

Note that the power actually contained within the angular area of 
a single K 3 pattern 

\E„(B h w )|W(l/«)Ml/e) 

is larger than the scattered power per mode by the factor 

50(l/e) 8<p(l/e) 
(X/6.) (X/6„) " 

This ratio gives a rough gaussian weighted measure of the number of 
modes that contribute to the intensity reaching a particular (0, <p). 

3.3.3. The spatial coherence function ot the slit-plane field 

In light-scattering experiments designed to extract spectral infor- 
mation from the scattered field using photocurrent correlation tech- 
niques, the feasibility of a particular measurement is critically de- 
pendent on the range of transverse spatial correlation that characterizes 
the observation plane field. 2 ' 3 - 27 The extent of the correlation is de- 



* See the discussion which follows eq. (113) and leads to eq. (124). 
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scribed quantitatively by the normalized mutual coherence function 39,40 

U ''' &,J {<|E.({ l i,;OI><|E.(f + A€ l , + A,; [«>}»' (100) 

where (£, 7;) and (£ + A£, 77 + A17) are two arbitrary points in the 
observation plane. The angular brackets denote an appropriate en- 
semble or time average. The function T(---) reaches its maximum 
value, T(- • ■) = 1, for A£ = A?7 = and, in general, decreases 
smoothly to zero as A£ and/or At? increase. The contour in A£ and A77 
around (£, 77) on which the coherence function reaches some specified 
numerical value may be taken as a measure of the area over which 
there is correlated temporal behavior of the two field amplitudes. 
When the main probe beam is derived from a source having perfect 
transverse spatial coherence, as is the case here, then the presence of 
the spatial incoherence in the scattered field is totally attributable to 
the scattering processes taking place in the illuminated volume. The 
spatial coherence properties of the scattered field are uniquely deter- 
mined at the exit face of the sample and are most easily specified 
analytically by calculating the mutual coherence function on a far- 
field reference sphere, 0, centered on the scattering volume. In purely 
formal terms, we can write 

r- / r o) - (E.(M)-E;(r + e,Q) 

7o(r ' e) -{(|E.(r,0| 2 )<|E.(r + e ,0| 2 )}^ (101) 

where r and r + p both terminate on the surface of the far-field sphere, 
0. Generally speaking, T (i, p) can be calculated in a straightforward 
fashion once it is assumed that the scattering perturbations satisfy 
certain basic stochastic criteria. 

The relationship between the observation plane coherence function 
r (£, V, A£, A77) and the far-field function T (r, p) depends, of course, 
on the detailed characteristics of the optical system which collects and 
images the scattered light, and must include the effects of aberrations. 
There are two alternative procedures that may be used to obtain this 
relationship. The first involves the use of the plane wave K, expansion 
of the scattering perturbations that was introduced in the beginning 
of this section. For the MK VI instrument, we have already calculated 
the slit-plane field produced by the scattering from the individual Ky. 
In the notation of eq. (89), we have 

E.a,T7;K y ) 

- E„ m(U „.) exp [ - S^L ] exp [ - -k^P ] , (102) 

where / is the effective focal length of the light collection system. In 
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theory, therefore, we could calculate T(£, V, A£, Ari) directly by ex- 
pressing the total slit-plane field as a sum over the E,(£, ij; Ky) and 
then performing the statistical average indicated in eq. (100). In the 
absence of aberrations or other imaging defects, this direct method 
represents the simplest approach. For perfect imaging, the individual 
plane-wave scattered field patterns are essentially nonoverlapping at 
the observation plane and the coherence function is effectively domi- 
nated by the contribution of a single Ky term. However, when imaging 
errors produce a significant overlap of the E„(£, 77; Ky) at the obser- 
vation plane, as is the case for the present apparatus, then obtaining 
the analytical form of T(£, 17; A£, Arj) by the direct method becomes 
a difficult mathematical problem. 

The second alternative approach involves a direct calculation of 
the far-field coherence function T (t, q) from which T(£, 17; A£, A17) is 
obtained by using the fundamental laws that govern the "propagation" 
of mutual coherence in an optical system. This latter method is 
generally the more useful when the light-collection system departs 
significantly from ideal imaging. 

For the scattering angles relevant to the MK VI instrument, the 
two-dimensional "phase sheet" model of the scattering sample may 
be used to simplify the calculation of T (r, q). For this two-dimensional 
model object, the reference sphere coherence function is given by the 
van Cittert-Zernike theorem 39 as 



, (103) 



To($, <p;e-d';<p- *>') 

2v j fdxdy\E (x, y) | 2 exp {ik [(d - d')x + {<p - <p')y~\] 

2tt J Jdxdy\E (x, y)\* 

where both r and p have been expressed in the cartesian angular 
coordinates d and <p. In eq. (103), the factor \E (x, y) | 2 is the illumi- 
nation function of the object, in our case the "phase-sheet" sample. 
The surface integral is to be taken over the entire (x, y) plane or over 
the open aperture of the object, as appropriate. It should be noted 
that the van Cittert-Zernike theorem will hold as long as the pertur- 
bations in the scattering "phase sheet" have a correlation distance, 
which is short compared to the characteristic spatial dimensions of 
\E (x, y)\ 2 . This condition is, in general, well satisfied in the typical 
scattering experiment. 

At small angles, where eq. (103) is valid, the far-field spatial co- 
herence function is independent of the absolute angular position of 
either observation point and depends only on the separations (0 — 6') 
and (<p — cp'). In terms of these difference variables, T (- • •) is just 
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the normalized Fourier transform of the source intensity. As such, it 
bears an extremely close resemblance to the instrumental profile 
calculated in Sections 3.1 and 3.2. 

For the MK VI instrument, the illumination function is the gaussian 

\E (x, y) \ 2 - El exp (- ^-±^ ) (104) 

and we have from eq. (103) 

/6x/2 
exp (-x 2 /a 2 ) exp [ik (d - d')x~]dx 

io W *;» - <r,v - *■; - -^ ^ 

/ exp ( — x 2 /a 2 )dx 

J-bx/2 
■by 12 



fb u l2 

/ exp ( — y 2 /a 2 ) exp \jik (<p — <p')y^dy 

J-b„/2 

/6„/2 
exp (-y 2 /a 2 )dy 
-i u /2 



X^^ j^ , (105) 



where b x and b v are the aperture dimensions at the scattering sample. 
As is evident from eq. (105), the coherence function factors for the 
case of gaussian illumination and we can write 

T (d, <p;6-e',<p- <p') = T' (Ad)T' (*<p), 

where Ad = 6 — 6' and A<p = <p — <p'. The functions T' Q (Ad) and 
T' (A<p) are given by the appropriate integrals in eq. (105). Each of 
these integrals is a finite domain Fourier transform of a gaussian kernel 
of the type considered in detail in Section 3.1 with respect to aper- 
ture apodization and vignetting. The only difference is that in eq. (105), 
the "intensity," exp ( — x 2 /a 2 ), replaces the "field," exp ( — x 2 /2a 2 ) f 
which appeared in the diffraction calculations. It is not hard to show 
that the factored coherence functions T' (A6) and To(A0) are identical 
to the normalized intensity profiles of Fig. 9 if one uses the corre- 
spondence 

IRWI-T^- ( 106 ) 

Given the form of the far-field reference sphere function T , we must 
now determine the relationship between T and the desired slit-plane 
correlation function. 

One of the fundamental results of coherence theory is that second- 
order mutual coherence functions, such as T( • • •), propagate according 
to the wave equations as "field" variables. That is, once T(- • •) is 
specified on any surface in an optical system, its form on any other 
surface in the system may be found by treating the coherence function 
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as one would any electric field distribution. Therefore, all of the usual 
wave-diffraction and/or geometrical-optics approaches used to analyze 
wave-front propagation in an optical system are directly applicable to 
the coherence function. 

For the MK VI apparatus, the coherence "field" described by T 
is identical to the far-field electric field that describes the instrument's 
directly transmitted probe beam, except for a numerical change in 
the beam-width parameter <r. The effective beam width which charac- 
terizes the mutual coherence "field," or, is related to the actual beam 
width of the instrument, a*, by the result 

(cr*) 2 



4 = *V-- ( 107 ) 



Except for this numerical change, the diffraction and aberration re- 
sults of Sections 3.1 and 3.2 may be used intact to describe the slit- 
plane coherence function. In terms of angular coordinates at the slit 
and the widths 50(l/e) and 8<p(l/e), which were used to characterize 
the instrumental profile, we have easily 

T(6, <p;d-d',<e- tp') = T'(M)r(* v ) 



= exp 



(6 - d'Y 
450* (1/e) 



exp 



(<P- <p') 2 
45*> 2 (l/e) 



(108) 

In the slit plane, as was the case on the surface of the far-field reference 
sphere, the slit-plane coherence functions are related to the intensity 
profile of the transmitted beam by the transformation 

Equation (108) is the basic result which may be used to evaluate the 
scattered or stray-light power-per-coherence region or estimate the 
number of coherence regions encompassed by a particular choice of 
main-slit size. For example, given the slit-plane scattered intensity 
\E.(d, <p)\ 2 , we can form the weighted integral 

dP j!?' y) = f fdd'd<p' \E.(tf, <p')\*LT'(6 - tf)T'{<p - <p')J, (109) 

which is a useful measure of the power-per-coherence solid angle as 
measured at the slit. 2 ' 327 In general, \E s (d', <p')\ 2 is slowly varying 
over the angular range where [T'(d - d')T'((p — </)] 2 is nonvanishing 
and can be removed from the integral to give 

dP ^ 6 ' v) = \E,(d, ^Mcoh^coh, (HO) 
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where the mean full-width coherence angles, A6coh and A<pcoh, are 
denned by the integrals 

&coh = &Qcoh&<pcoh 

= f[T'(A6)Jd{Ad) flT'(A<p)yd(A<p). (Ill) 

Combining eqs. (108) and (109) gives for the MK VI instrument 

aIcoh = \&0(l/e) = 24.1 M rad 
A~vcoh = V27r5<p(l/e) = 140.4 /xrad 

and eq. (110) becomes 

dP ;^ v) = \E.(e, *)|W2rt*(l/e)V2SMl/e). (113) 

Earlier in this section, we obtained an expression for | E,(d, <p) | 2 based 
on a plane-wave-mode expansion of the scattering perturbations. That 
result may be used in eq. (113) to yield a relationship between the 
observed scattered power-per-coherence solid angle and the scattered 
power-per-Ky mode. From eqs. (99) and (113), we find 

dP,(d h ip,) _ V2^fl(l/e) Vgjgg(l/g) p ,„, m4 v 

dtlcoH ~ (X/6.) X (X/K) ^* ih U14; 

The product of the correction factors 



l2vM(l/e) V27r5<p(l/e) 



(X/6.) ~ (\/b y ) 

is a rough measure of the number of modes that contribute to the 
power observed in a single coherence region at the slit plane, while the 
individual terms indicate the extent of the multiple mode contribution 
in the 6 and directions. For the MK VI instrument at full aperture, 
the numerical values of the correction factors are 



V2~^ 5fl(l/e) 
(X/bx) 

V2^r Ml/g) 
(X/6y) 



= 2.4 

(115) 
= 14.0. 



The results given in eqs. (99), (113), and (114) together with the 
known form of the instrumental profile may be combined in various 
ways to calculate normalized scattering cross sections from measured 
slit-plane intensities. One important calculation of this type is to 
express the observed stray-light levels in the MK VI apparatus in 
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terms of an equivalent scattering cross section. For the experimental 
profile curves shown in Fig. 4, we may write down an analytical ex- 
pression for the measured stray-light photocurrent, i(d), as 

/ \E, t (0, <p)\'d0d Vt (116) 

-AB BL /2 J-A faL /2 

where \E tl {6, <p)\ 2 is the stray-light intensity at the slit plane and 
AOsl and A<p S L specify the full-width slit dimensions in angular units. 
The proportionality factor a relates the photocurrent to the optical 
power passed by the slit and includes the detector quantum efficiency, 
light-collection losses, etc. If the intensity \E, t (d, <p)\ 2 is relatively 
constant over the slit aperture, we have simply 

i(d) = a\E st {6, 0) \*A6sl&<Psl, (117) 

where we have assumed that A<psl is situated symmetrically around 
<p = 0. Combining this result with eq. (110) gives the relation between 
the measured photocurrent and the stray-light power-per-coherence 
solid angle as 

m -a*g&g*pB-. (118) 

aUcoH AOcoh A<pcoh 

To eliminate the unknown proportionality constant a, we make use of 
photocurrent observed at 6 = 0, the peak of the directly transmitted 
beam. Given the normalized slit-plane intensity profile of the direct 
beam, 1 (0, <p)//(0, 0), we can calculate the fraction of the total beam 
power, P , passed by the slit at 6 = as 

■A9S1/2 [*vslI2 /(fl, g) 

-AySL/2 - f (0, 0) _ (119) 



-A9SL/2 J- 



If. 



mi™* 

li a I (0, 0) 



where for the MK VI apparatus we have 

ULA - eX n T- -±- I exp I"- * 1- 

m o) - exp L mv*) J p L wo./*) J 

The numerical value of the error function integrals in eq. (119) could 
be obtained from tabulated results for particular values of A6 S l and 
A(psl\ however, in the present case where the slit dimensions satisfy 
the inequalities 

A6sl « 80(l/e) (l2Q) 

A^ 5L »5^(l/e), k 
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we have the more useful analytical result 



Ad SL 

y -^se(i/e)' (121) 

The measured peak photocurrent, i(0), is then 

Dividing eq. (118) by eq. (122) gives the useful result 

i{6) = 1 dP at (d,0) 1 A<p SL 

t(0) Po rfficoH ^A^pcoh' ( ' 

If desired, the quantity dP, t {d, 0)/dti C oH can be replaced with the 
stray-light power per mode, P„(Kj), by using eq. (114). This gives the 
very useful relationship 

i(0) (P.,(K,-) V&r"gfl(l/e) A<p SL 

t(0) " (Po (X/6.) (\/6„)' KLZ * } 

IV. EMPIRICAL OBSERVATIONS ON THE STRAY-LIGHT BEHAVIOR 
OF OPTICAL ELEMENTS AT VERY SMALL ANGLES 

Very little information of a quantitative nature is available con- 
cerning the imperfection scattering of optical elements at very small 
angles. As a result, the design and testing process leading to the 
present vsa instrument involved a significant amount of trial and error 
evaluation of various optical systems in a search for the desired stray- 
light performance. During this process, a certain amount of empirical 
information was obtained relating to the imperfection-scattering ques- 
tion. This section presents a brief discussion of these observations and 
their influence on the configuration adopted for the MK VI instrument. 

4.1 Reflecting versus refracting optics 

It is clear from a comparison of Figs. 1 and 6 that the implementation 
of a vsa scattering instrument using lenses would be significantly less 
involved than the A IK VI off -axis mirror arrangement. The refracting 
system also has the advantage of strictly zero off-axis aberrations 
(coma, astigmatism, and distortion), although a "best form" single- 
element lens does have eight times the spherical aberration of an 
equivalent spherical mirror. 37 In fact, the earliest version of the present 
apparatus utilized precisely the kind of "straight-through" lens system 
illustrated in Fig. 6. This arrangement was abandoned because of 
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two problems : 

(i) The presence of Newton's interference fringes crossing the 

illuminated field. 
(ii) An excessive stray-light background. 

The first problem arises because of the partial reflectivity of the two 
lens surfaces and can be solved to some extent through the use of anti- 
reflection (ar) coatings. However, even the best antireflection coated 
lens will form far-field Newton fringes with an integrated inten- 
sity of about j percent of the incident beam power. This fact makes 
the refracting components generally unacceptable in a vsa system. 
The presence of these extraneous reflections and their associated 
interference fringes creates an intense fixed-pattern nuisance back- 
ground which can make it impossible to observe the angular de- 
pendence of the sample scattered light. The stray-light background 
problem is a manifestation of small-angle scattering at the lens which 
may originate from three possible sources : 

(t) Lens surface "roughness" or nonconformity (at least two 

surfaces) . 
(ii) Index of refraction inhomogeneity in the lens bulk material. 
(Hi) ar coating thickness nonuniformity (at least two surfaces). 

By way of comparison, the possible sources of imperfection scattering 
from a first-surface reflector are 

(i) Mirror surface "roughness" or nonconformity (one surface). 
(it) Reflective coating (s) thickness nonuniformity. 
(Hi) Reflective coating reflectivity nonuniformity. 

From a theoretical standpoint, one should be able to evaluate the 
seriousness of each of these defects a priori by calculating the surface 
and/or bulk inhomogeneity scattering. This calculation is straight- 
forward if one has available the spatial form of the roughness in terms 
of the spatial correlation function and the rms roughness amplitude. 
The effect of roughness or inhomogeneity is to impose a spatially 
random-phase perturbation in the optical path. The scattering that 
takes place as a result of this perturbation can be calculated via the 
same "phase-object" approach which is used for the primary scattering 
sample (see Section 2.6). The stray-light intensity observed at some 
specified scattering angles 6 and <p is given by the Fourier transform 
of the roughness correlation function at a wave vector |K| = 2t/A 
satisfying the appropriate small-angle kinematic conditions. Un- 
fortunately, the roughness wavelengths corresponding to the angular 
range of interest here (10 -3 cm < A < 1 cm) are determined by a 
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spatial region of the roughness correlation function about which very 
little is presently known. This wavelength regime presents difficult 
measurement problems and is generally not probed by conventional 
roughness-testing techniques. The data that is available comes from 
two measurement techniques that tend to flank this regime on the 
short and long wavelength sides : 

(i) The feco* interferometer and allied methods 41 - 42 that exhibit 
good surface deviation resolution, 1 A to 10 A, but are useful 
only at short wavelengths (A < 1000 A). 
(it) Conventional "surface-conformity" techniques such as the 
Foucault knife-edge and Twyman-Green interferometer tests 
that are useful primarily at longer roughness wavelengths 
(0.1 cm to 100 cm) and which exhibit relatively poor surface 
deviation resolution (50 A — » 2000 A). 

The stray-light measurements that were made during the course of 
the evolution of the present instrument provided the most sensitive 
roughness and inhomogeneity test for this awkward wavelength range. 
It was found experimentally that, for lenses and mirrors of the same 
fraction of the "state-of-the-art," the stray-light level of a refracting 
instrument was roughly 20 times that of its reflecting counterpart. In 
neither case did the vsa stray-light level correlate well with known 
short wavelength roughness and inhomogeneity data. Both types of 
components exhibited a spatial roughness spectrum that was strongly 
enhanced at long wavelengths. This enhancement did not appear to 
depend as strongly on the "surface-figure" of the component as one 
might be led to expect by qualitative theoretical arguments. 

Comparisons were also made between mirror components having 
multilayer dielectric coatings and those with a conventional SiO- 
protected aluminized surface. The aluminized coatings can suffer from 
a spatially varying reflectivity caused by surface oxidation while high- 
reflectivity dielectric films tend to have a significantly smaller reflec- 
tivity modulation. However, the stray-light measurements showed no 
significant difference between the two types of coatings on similar 
"quality" substrates. Apparently the cumulative roughness of the 
greater number of dielectric layers offsets the dielectric coating's 
potential advantage. 

4.2 Main scanning silt 

Another major contributor to the stray-light level in earlier versions 
of the MK VI apparatus was the main angle-scanning slit. The slit 
selected for this application is a commercial Spex unit normally used 

* Fringes of Equal Chromatic Order. 
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as an intermediate or exit slit on a double-grating spectrometer known 
for its low stray-light background. This fact not withstanding, severe 
stray-light problems were encountered in predecessors of the MK VI 
that had this slit located directly at the focal plane of the collecting 
mirror (mirror M 6 shown in Fig. 1). The origin of this problem was 
traced to scattering of the direct beam by the beveled surfaces of the 
slit jaws and to quasi-specular reflection from the slightly flattened 
and rounded jaw edges. This source of background by itself was of 
sufficient intensity to completely swamp the sum total of all other 
stray-light sources in the instrument. 

This problem was solved in the MK VI apparatus by occulting the 
directly transmitted beam, before it reached the main scanning slit, 
with a precision knife-edge fabricated of highly attenuating black glass 
plate. The use of glass instead of metal permits the edge defining 
surfaces to be optically polished without cold flaw and rounding. In 
addition, the included angle formed by the edge surfaces is made 
obtuse, rather than the acute angle normally used, to avoid the 
feathering problems and surface irregularity enhancement associated 
with small included angles. The salient geometrical features of the 
knife-edge are illustrated in Fig. 20. The actual occulting edge is formed 
by a single beveling operation on polished flat stock and is oriented in 
use such that an incoming ray strikes the beveled face at the quasi- 
Brewster angle. The beveling angle is chosen so that the ray which is 
refracted into the plate travels parallel to the plate surfaces and is 
totally absorbed. 

The improvement in stray-light level obtained by using the knife-edge 
to occult the direct beam, rather than relying solely on the main scan- 
ning slit, can be seen in Figs. 3 and 4. The improvement amounts to 
roughly an order of magnitude over the angular range of interest. 

4.3 Aberration corrections and stray light 

The reader familiar with optical system design will realize that the 
aberrations present in the MK VI instrument could be "corrected" 
using well-known techniques. However, the application of these cor- 
rection methods has two drawbacks: cost and reduced stray-light 
performance. The simplest corrective measures, those which add the 
fewest number of optical elements to the basic apparatus, entail the 
use of off-axis fabricated, aspheric reflecting and/or refracting elements. 
These types of elements are, in general, exceedingly costly to fabricate. 
More sophisticated aberration-corrective designs, utilizing only spheri- 
cal optics, require a larger number of additional elements. In either 
case, of course, the presence of additional optical surfaces means 
degraded stray-light performance. Furthermore, any corrective design 
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Fig. 20 — Geometrical features of the occulting knife-edge used in the MK VI 
apparatus. 

relying on the use of refracting elements will be further penalized by 
the excessive small-angle stray light which these elements generate. 

APPENDIX A 

Finite Slit-Width Effects in the Scanning of Gaussian Intensity Profiles 

When a gaussian focal-plane profile is scanned by a finite-width slit, 
the transmitted power is proportional to the integral 

/(fe A) = / exp [- e/5e(l/e)ldi-, (125) 

where £ gives the position of the center of a slit whose width is 2A. 
By writing the spatial coordinate £ as 

£ = £o + r (126) 

and making a change of variable, we may put eq. (125) into the form 
J (to, A) = exp[-#/8?(l/e)] 

X J* exp [- (2£„ + r 2 )/5^(l/e)]df . (127) 
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Table IX — Angular displacement, d , at which an unbroadened 

gaussian and a slit-broadened gaussian reach specified 

fractions of peak intensity 
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38.678 


0.36 
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For reasonably small values of the ratio [A/5£(l/e)], the gaussian term 
in the integrand of eq. (127) may be approximated by the leading 
term in its Taylor's series expansion 



exp[-r 2 /8^(l/e)] = 1 - 



WOT*) 



+ 



with a maximum error exp [— A 2 /5^(l/e)]. Within this approximation, 
the remaining integral can be calculated in a straightforward manner 
to give 



«/(£„, A) = (2A) exp [-£o725£ 2 (lA)] 



sinh [2£ A/5g 2 (l/e)] 
[2€oA/«?(l/e)] 



(128) 



Since the function (sinh x)/x tends to unity as x goes to zero, the 
normalized slit-broadened profile is 



7(07A) " GXP L &/ ** (1/6)J [2€oA/Wl/<0] 



(129) 



or its equivalent written in terms of the scattering angle 6 = £//. 
Clearly, in the limit A — ► 0, eq. (129) describes the correct unbroadened 
gaussian. For A^O, the principal effect of the (sinh x)/x correction 
term is to push up the tails of the profile while leaving the peak of the 
gaussian relatively unaffected. A good quantitative feeling for the 
nature of this correction may be obtained by solving for the off-zero 
displacements, £0, at which the broadened and unbroadened profiles 
reach specified fractions of their peak intensity. These £0 values then 
specify the profile half-widths at the corresponding intensity level. 
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For the curves presented in Section II, the relevant numerical 
parameters, expressed in angular units, are : 

80(1 /e) = (l//)5£(l/e) = 9.6 M rad 
A fl = (1//)A = 1.0 M rad. 

Table IX gives the calculated half- width values O = (W/) obtained 
from eq. (129) for various choices of the ratio J(6 , A)/J(0, A). For 
comparison, the table also lists the corresponding half-widths of the 
uncorrected gaussian, and the percentage of line-width increase caused 
by the slit-width correction. As is evident from these results, the 
effect of the (1/a:) sinh x correction term is to alter the gaussian profile 
in such a way that the observed half-widths are an essentially constant 
■percentage larger than the true values. 

APPENDIX B 

Numerical Evaluation of the Diffraction Profile of Apertured 
Gaussian Illumination 

Equation (47) gives the basic integral for the truncated gaussian 
diffraction profile as 

J?„ rft/2 

E W " 77TTT / ex P (-*/•") ex P P(VAo)*. (130) 
UMP J -bit 

This expression may be put into a form more suited to numerical 
computation as follows. Wo write the exp «'(• • •) term as 

exp [i{2ir/ f\ Q )$,x~\ = cos Kx + i sin Kx 
with 

K = 2rf _ 2*6 (131) 

JAo Ao 

and note that the sin Kx integral vanishes by symmetry. Next by a 
change of variable 

* = 6 ^, (132) 



we obtain 



E(i) = -M^ r* 1 cos cw exp (-aho 2 )dw, (133) 

(JAo)' J w=0 



where a and c are denned as 
Kb 






= {f^)2 = M> a= S?' (134) 

The gaussian in the integrand is now expressed in terms of its Taylor's 
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series expansion 

exp (-a?x 2 ) = D -j 

n=0 W! 

to give E(£) as 

fl(0 = Jg* £ ( ~ 1)n . (a2) " f w«»co8 (cw) dw. (135) 

Equation (135) forms the basis for the numerical computation of the 
profiles. 

A simple closed-function form for the w integrals in eq. (135) does 
not exist; however, recursive relations among these integrals can be 
found from the standard integrals 

/ x m cos ex dx = / x m ~ l sin ex dx (136) 

Jo c c Jo 

and 

f 1 _ i • i cose . (m + 1) f 1 __- , , 1Q7 \ 

/ x m ~ l sin ex da; = \- - / x m ~ 2 cos ex dx. (137) 

Jo c C Jo 

Defining 

L m (c) = I x m cos ex dx, (138) 

./o 

we easily obtain the following recursion formulae from eqs. (136) 
and (137) : 

/ . oxr ,n /m + 3\ . . (ro + 3)(ro + 2) 

(m + 3)Lm+a(c) = ( ! — J sin c + ^ ^ - cos c 

(m + 3)(m + 2) r , , 1NT , >.-, / 1Qft \ 
- ^ ^ [(w + l)L«(c) j (139) 

(m — l)L m _ 2 (c) = f — J sin c + 



' )S c 



m(m + 1) 
From eq. (136), we also have for m = 



[(m + l)L m (c)]. (140) 



L (c) = S££. (141) 

In terms of the L m (c), the expression for the diffracted field takes the 
series form 

JS(© = ^ f ( - 1)n(a T^" (C) d42) 

U A 0.r n = W! 
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and the desired normalized intensity profile is 

/(0) 



f - (-l)«(a 2 )»L 2n (c) 


2 


1 n% «! 




f (-l)»(a 2 )»L 2n (0) 


2 


nt-o ft! 





(143) 



The zero argument L m 's can be written down explicitly from eq. 
(138), viz. 

MO) = ztT (144) 

For the numerical results reported here, the series in eq. (143) were 
truncated at some n = Umax by testing the value of (l/n!)(a 2 ) n L 2n (0) 
and terminating when this quantity was smaller than some chosen 
convergence criterion, e. In the present case, e was set at e = 10 -12 . 
For the largest (b/a) value, (b/a) = 8.33, where the gaussian kernel 
of eq. (133) is 

exp (-oV) = exp (-8.68m; 2 ), 

49 terms in the series were required for convergence. 

For each individual pair of values for c and Umax, the required string 
of L m 's are generated by two subroutine programs. 

B.1 Subroutine No. 1, c < 1 

When the quantity c = (Kb/2) is less than one, the L m 's are obtained 
by the following procedure. 

(i) Calculate L m (c) for m = 2n M Ax directly from the defining 
equation (138), using the Taylor expansion for cos ex to write 



» /_lW2n r\ 



_ « (-DV a45) 

" n tr (2tt)!(m + 2n+l) k ' 

(ii) Truncate the sum in eq. (145) when c 2n /(2n) \(m + 2n + 1) is 

less than 10~ 12 . 
(m) Use this result for L n u AX (c) to obtain the required L m 's via 

the backward recursion formula, eq. (140). 

B.2 Subroutine No. 2, c> 1 

When the quantity c = (Kb/2) is greater than 1, the L m 's are found 
by a two-part procedure that depends on the value of timax. 

(i) For m = 2n values for which the inequality m = 2n < c is 
satisfied, use L\(c) = (sin c)/c and the forward recursion rela- 
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tion, eq. (140). If 2wmax is less than c this first step gives 
all required L m 's. 
(ii) If 2?i M ax is greater than c, set 

r , x _ . cos c 
L m (c) 



r^j 



m + 1 



for some m ~2> 27i MA x and work backward using recursion relation eq. 
(140). The calculated string of L m 's is joined onto the forward recursion 
values from step (i) for some m tt c and then renormalized. 

This rather elaborate procedure for calculating the L m 's is made 
necessary by the rapid accumulation of numerical round-off errors 
which arise in the repetitive application of the basic recursion formulae. 
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